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FOREWORD 

As  part  of  the  on-golng  program  in  "Decision  and  Control 
Models  in  Operations  Research,"  Mr.  Richard  Ehrhardt  has  developed 
a number  of  formulae  that  give  approximate  values  for  the  operating 
characteristics  of  several  classes  of  (s , S) / Inventory  policies. 

The  approximations  are  accurate  and  require  substantially  less 
ct)mputatlona 1 effort  than  exact  calculation  of  the  operating 
characteristics.  Many  of  the  approximations  are  in  a form  that 
make  them  particularly  amenable  to  sensitivity  analysis  and  systems 
design  studies. 

I Other  related  reports  dealing  with  this  research  program  are 

given  below. 

Harvey  M,  Wagner 
Principal  Investigator 
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ABSTRACT 


OPKRATING  CHARACTERISTIC  APPROXIMATIONS  FOR  THE 
ANALYSIS  OF  (s,S)  INVENTORY  SYSTEMS 

Richard  Ehrhardt 

University  North  Carolina  - 1977 

The  operating  characteristics  of  (s,S)  inventory  systems  are  often 
difficult  to  compute,  making  sensitivity  analysis  a tedious  and  often 
expensive  undertaking.  Approximate  expressions  for  operating  character- 
istics are  presented  with  a view  towards  simplified  sensitivity  analysis. 

The  operating  characteristics  under  consideration  are  the  expected 
values  of;  Total  cost  per  period,  period-end  Inventory,  period-end 
stockout  quantity,  replenishment  cost  per  period,  and  backlog  frequency. 
The  approximations  are  obtained  by  using  least-squares  regression  to  fit 
simple  functions  to  the  operating  characteristics  of  a large  number  of 
parameter  s<.-t tings.  Accuracy  to  within  a few  percent  of  actual  values 
is  typical  for  most  of  the  approximations. 

Potential  uses  of  tne  approximations  are  illustrated  for  several 
idealized  design  problems. 
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1.  (s,S)  INVENTORY  POLICY  OPERATING  CHARACTERISTICS 


There  are  many  situations  In  which  an  Inventory  systems  designer  can 
use  estimates  of  operating  characteristics  of  the  system.  For  example, 
management  may  desire  forecasts  of  invetitory-on-hand , or  system  operating 
costs.  The  systems  designer  may  also  confront  questions  such  as 

(a)  Wliat  are  the  cost  and  service  implications  of  consolidating 
demand  from  several  distinct  warehouses  into  a single  central 
warehouse? 

(b)  How  do  parameter  changes  affect  the  system's  operating 
characteristics? 

(c)  By  how  much  do  costs  rise  when  service  is  increased? 

(d)  What  are  the  effects  of  changing  the  review-period  length? 

To  facilitate  analysis  of  design  issues  of  this  type,  we  seek,  simple 
approximations  for  the  following  operating  characteristics:  average  hold- 
ing cost  per  period,  average  backlog  cost  per  period,  frequency  of  periods 
without  backlogs,  average  replenishment  cost  per  period,  and  average  total 
cost  per  period.  These  characteristics  are  defined  mathematically  in 
Section  2.1  of  this  report. 

1 . 1 The  Model 

Tliroughout  tills  report,  we  deal  with  a single-item  inventory  model. 

We  assume  periodic  review  of  an  item's  Inventory  level  and  employ  a station- 
ary, discrete-time  stocliastlc  process  to  describe  the  Item's  demand.  The 
demand  sequence  consists  of  independent,  identically  distri- 

buted random  variables  taking  non-negative  integer  values. 
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Demands  are  met  as  long  as  stock  on  hand  is  sufficient;  when  a stock-  * 

out  occurs,  the  unfilled  demand  is  completely  backlogged  until  a stock  i 


replenishment  eventually  arrives. 


Items  kept  In  Inventory  are  assumed  to  be  conserved,  there  being  no 
losses  by  deterioration,  obsolescence,  or  pilferage;  disposal  Is  not 
allowed.  Invenlory-on-hand  at  the  end  of  a period  Is  the  inventory  from 
the  previous  period  plus  any  replenishment  that  arrives,  less  demand. 
Negative  inventory-on-hand  represents  the  amount  of  backlogged  demand. 
Replenishments  are  assumed  to  be  delivered  a fixed  lead  time  L periods 
after  being  ordered.  The  time  sequence  of  events  In  any  period  is  taken 
to  be  order,  delivery,  demand. 

We  assume?  no  time  discounting  of  costs  and  postulate  an  unbounded 
horizon  over  which  the  item  is  demanded  and  stocked.  We  seek  to  minimize 
expected  total  cost  per  period. 

The  cost  of  a replenisliment  quantity  q is  assumed  linear  with  fixed 
ordering  cost  K and  constant  unit  cost  c 


c(q) 


|k  + cq 

I 0 


for  q > 0 
for  q = 0 . 


Items  are  not  lost  from  Inventory,  so  demand  is  completely  filled. 
Since  costs  are  not  discounted,  the  constant  unit  cost  c is  not  a factor 
in  choosing  a minimum  cost  policy,  and  is  suppressed  hereafter. 

The  inventory  holding  cost  is  proportional  to  any  stock  on  hand  at 
unit  cost  h 

f hi  for  i > 0 

h(i)  = 

lo  for  i •=  0 , 

and  the  unit  penalty  cost  p is  applied  to  any  quantity  on  backorder 
at  the  end  of  each  period 
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for  i > 0 
for  i < 0 . 

The  resulting  total  cost  function,  therefore,  is  linear  in  K , p , 
and  h and  wo  may  scale  these  parameters  so  that  the  value  of  the  unit 
holding  cost  h is  unity.  Non-trlvial  changes  in  cost  arise  only  with 
ciianges  in  the  ratios  K/h  and  p/h  . 

1 . 2 Inventory  Replenishment  Policies 

We  postulate  that  control  over  replenishment  Is  exercised  by  an  (s,S) 
policy:  whenever  inventory  x on  hand  and  on  order  at  the  start  of  a 
period  drops  below  the  value  "little  s,"  an  order  is  placed  for  a replenish 
merit  of  size  S-x. 

Given  our  assumptions,  when  the  demand  distribution  and  the  economic 
parameters  are  known,  there  is  an  optimal  policy  that  has  the  (s,S)  form 
[iglehart  (1963a, b),  Veinott  & Wagner  (1965)].  When  tlie  demand  distribu- 
tion is  not  known,  even  thougli  this  is  the  only  assumption  relaxed,  an  optl 
raal  policy  may  no  longer  be  of  the  (s,S)  form.  Nevertheless,  in  this  ex- 
periment we  employ  an  (s,S)  policy,  since  it  is  in  popular  use  in  the 
applied  situation  of  incomplete  Information. 

We  seek  approximations  for  the  operating  characteristics  of  three 
(s,S)  policy  rules:  optimal  policies.  Power  Approximation  policies,  and 
Statistical  Power  Approximation  policies.  For  each  policy  rule  we  develop 
simplified  functions  for  approximating  the  operating  characteristics.  We 
then  fit  the  parameter  values  of  the  functions  to  the  observed  cliar.ic  ter  Is- 
tlcs  of  a multi-item  system  using  least-squares  regression. 


1.2.1  Opt  iimi  1 Pol  Icles 

We  develop  approximate  expressions  for  optimal  policy  operating 
cltarac t or  1 St  1 cs  in  Sections  2 and  3.  The  approximations  are  compared 
with  actual  values  of  the  operating  characteristics,  calculated  with 
the  algorithm  of  Veinott  and  Wagner  (1965). 

1.2.2  The  Power  Approxlmat ion 

The  Power  Approximation  [Ehrhardt  (1976)]  is  an  algorithm  for  com- 
puting approximately  optimal  values  for  (s,S)  using  only  the  mean  p 
2 

and  variance  a of  demand.  The  algorithm  is  executed  as  follows. 

Let 

(1)  D = (1.463)M'^^'^(K/h)'^^®[(L+l)o^]'°^^^ 

p 

and 

Sj  = (L41)  p + [(L+l)p]'^^^  (o^/u)’^^^  U(z) 

(2) 

^1  ' “l  ^ ' 

where  U(z)  is  given  by 


U(z)  » . 182/z  + 1.142  - 3.466Z  , 

. . J iLi!V/hr^ 1 '' 

1 (l-H^)  [(L+l)o'’]-^^^l 


If  D /p 
P 


is  greater  than  1.5,  let 


Otherwise 


compute 

(4)  - (L+1)  p + V [(L+l)o^]**  , 
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whert-  V Is  Lhi>  solution  to 


('’)  <t’(v)  = p/(p+h), 

and  4'(')  is  the  cumulative  distribution  function  of  the  unit  normal 
distribution.  The  policy  parameters  are  then  given  by 
s = minimum  {Sj  , S^} 

(6) 

S = minimum  {S^  , S^}  . 

if  demands  are  integer  valued,  , and  S2  are  rounded  to  tiie 

nearest  integer. 

Approximate  expressions  for  Power  Approximation  operating  cliaracter- 
istics  are  developed  in  Sections  2 and  3.  The  approximations  are  compared 
with  actual  values  of  the  characteristics,  calculated  using  the  methods  of 
Veinott  and  Wagner  (1965). 

1.2.3  The  Stalls tical  Power  Approximation 

Of  course,  in  real  applications,  tlie  mean  and  variance  of  demand  are 
not  always  known.  For  the  situation  where  only  sample  statistics  of  pre- 
vious demands  are  available,  tlie  Statistical  Power  Approximation  described 
below  can  be  implemented.  This  decision  rule  implicitly  assumes  that  de- 
mand is  stationary.  Such  a rule  derived  for  stationary  conditions  may  be 
a leasonable  approximation  to  an  optimal  rule  when  the  demand  process  is 
mildly  non-stat ionary,  provided  that  the  policy  parameters  s and  S are 
revised  periodically  to  meet  the  changing  conditions. 

We  assume  in  this  study  that  a denujnd  history  of  fixed  length  is  kept 
to  make  each  revision,  and  equal  weight  Is  given  to  each  observation.  This 
is  not  optimal  If  the  demand  process  is  known  to  be  stationary,  for  then 
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the  entire  history  should  be  accumulated  to  give  progressively  better 
knowledge  and  performance.  Even  when  denuind  is  known  to  be  non-station- 
ary,  but  varying  in  a regular  manner,  such  as  by  a trend  or  periodic  cycle, 
or  both,  an  optimal  decision  rule  would  generally  utilize  the  entire 
h i story . 

The  decision-maker  usually  is  not  in  a position  to  know,  however, 
that  conditions  observed,  even  over  the  entire  past  history,  will  continue 
to  prevail.  This  provides  justification  for  making  frequent  revisions, 
placing  greater  weight  on  observations  from  the  immediate  past  and  less  on 
earlier  history.  For  tliis  study  the  admittedly  arbitrary  choice  has  been 
iiuide  to  keep  a history  of  fixed  length  and  give  equal  weight  to  all  obser- 
vations in  this  history.  Let  T be  the  number  of  periods  between  policy 
revisions,  which  will  be  termed  the  revision  interval;  assume  that  a lils- 
tory  of  T periods'  demands  is  kept  for  use  at  each  revision. 

The  statistics  required  by  our  decision  rules  are  the  sample  mean 
and  variance  of  demand,  C and  v , respectively.  If  t is  a period  at 
the  beginning  of  which  revision  is  made,  then 

Wlien  using  the  Statistical  Power  Approximation  we  periodically  obtain 

_ _ 2 

values  for  (s,S)  by  substituting  C and  v for  p and  o in 
equations  (1)  through  (6). 

We  develop  approximate  expressions  for  Statistical  Power  Approxi- 
mat ion  operating  characteristics  in  Section  3.  The  approximations  are 
compared  with  cstintites  of  actual  values  of  the  characteristics,  obtained 


from  computer  simulation  experiments. 
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1.3  Methods  of  Approximation 


In  Section  2 we  develop  "analytic"  approximations  for  the  operating 
cliaracterlstics  of  optimal  and  Power  Approximation  policies.  We 
approximate  the  exact  expressions  with  simplified  functions,  generalize 
these  functions,  and  then  fit  their  parameter  values  to  the  observed 
characteristics  of  576  policies  using  least-squares  regression.  The  re- 
sulting approximations  tend  to  be  extremely  accurate  but  rather  diffi- 

I 

cult  to  evaluate. 

In  Section  3 we  develop  "multiplicative"  approximations  for  fixed 
and  statistical  policies.  The  approximations  are  based  on  simple  function- 
al forms  that  are  motivated  by  empirical  observations.  The  functions  are 
multipl Icatlvely  separable  into  simple  expressions  that  each  depend  on 
two  variables  at  most.  We  use  least-squares  regression  to  fit  the  para- 
meters of  these  functions  to  the  observed  characteristics  of  576  fixed 
policies  and  288  statistical  policies.  The  resulting  approximations  are 
accurate  and  easy  to  evaluate. 

In  Section  4 we  use  the  approximations  to  study  sensitivity  analysis 
questions  that  typically  confront  systems  designers.  We  also  use  the 
approximations  to  draw  graphs  displaying  the  operating  characteristics  as 
functions  of  system  parameter  settings. 
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2.  ANALYTIC  APPROXIMATIONS  FOR  FIXED  POLICIES 


In  this  section  we  develop  approximate  expressions  for  the  operating 
characteristics  of  fixed  (specified)  policies.  We  approximate  exact  ana- 
lytic expressions  with  simplified  functions,  generalize  these  functions, 
and  then  fit  their  parameter  values  to  the  observed  characteristics  of 
576  items  using  least-squares  regression.  The  576-ltem  system  is  formed 
by  using  a full-factorial  combination  of  the  parameters  in  Table  2.1. 


Table  2.1 
System  Parameters 


Factor 

Levels 

Number 
of  Levels 

Demand  distribution 

Poisson  (o^/y  " 1) 

Negative  Binomial  (a^/p  “ 3) 
Negative  Binomial  (o^/y  “ 9) 

3 

Mean  demand 

2,  4,  8,  16 

4 

Replenishment  leadtime 

0,  2,  4 

3 

Replenishment  setup  cost 

32.  64 

2 

Unit  penalty  cost 

4,  9,  24,  99 

4 

Unit  holding  cost 

1 

1 

Policy 

Optimal  policy. 

Power  approximation  policy 

2 

The  approximations  in  this  section  are  functions  of  only  tlie  economic 
parameters,  policy  parameters,  and  the  mean  and  variance  of  demand.  We 
obtain  results  for  holding  cost,  backlog  cost,  backlog  protection,  replenish- 
ment cost,  and  total  cost.  On  the  average,  the  approximations  deviate  from 
the  actual  values  of  these  characteristics  by  0.7%,  A.1%,  0.7%,  0.1%,  and 
1.9%  respectively. 
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2 . 1 The  Analytic  Derivation 


Consider  the  model  of  Section  1.1  and  assume  that  demand  follows  a 
probability  density  ((iC*)  and  cumulative  distribution  4>(‘).  Let 
and  $ (•)  be  the  n-fold  convolutions  of  these  functions.  We  consider 

the  following  operating  characteristics  of  fixed,  infinite-horizon  (s,S) 
policies 

H = average  holding  cost  per  period 
B = average  backlog  cost  per  period 
(7)  P E backlog  protection,  i.e.,  frequency  of  periods  without 

backlogs 

R E average  replenishment  cost  per  period 
T E average  total  cost  per  period. 

Let 

muu  n 

= ^n=l 

n=l  ^ • 

We  have,  as  in  Roberts  (1962),  the  exact  relationships 

H » h[l-Hi(D)]"^{/^/^”^(S-y-x)(ti*^^‘'’^^  (x)m(y)dxdy 

+ /^(S-x)((i*^^^^^  (x)dx) 

' o 

B = p[H/h  + (L+l)u-S]  + p[l+M(D)]"^/°ym(y)dy 
(8)  P • (l+M(D)]"^/V^^'*'^^(S-y)m(y)dy  + 4>* (S)  } 

R - K[1+M(D)]"^ 

T - H + B + R , 


A 
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t 

where 

D : S - s . 

It  Is  difficult  to  obtain  any  insights  from  (8)  regarding  the  sensitivity 
of  the  operating  characteristics  to  values  of  model  parameters.  indeed 
it  is  exceedingly  complicated  just  to  calculate  values  of  the  characteris- 
tics for  a given  set  of  parameter  values.  We  proceed  to  simplify  the 
lorm  of  expressions  (8)  by  Introducing  approximations  for  the  functions 
in(*)  , M(*)  , and  iji  . 

Replenishment  frequency  In  (8)  is  given  by  [l+M(D)]  To 
approximate  M(*)  we  use  the  following  result  of  Smith  (1954). 

M(x)  = x/p  + o2/(2g2)  - 1/2  + o(l)  , x^®  . 

1 

i 

This  yields  the  approximate  value  for  replenishment  frequency 

(9)  [1+M(D)]'^  i g/[D  + (p4o2/ij)/2]  H p . 

We  identify  the  quantity  (S-y)  in  (8)  as  inventory-on-hand-plus-on- 
order  (after  ordering),  with  stationary  distribution  F(*)  given  by 

I M(y)/[1+M(D)]  , s ^ S-y<S 

(10)  F(S-y)  - 

( 1 , S-y  = S . 

The  probability  density  f(‘)  of  inventory-on-hand  plus-on-order  (after 
ordering)  on  the  open  Interval  [s.S)  is 

f(S-y)  ■ m(y)/ (1+M(D) ] . 

We  approximate  f(*)  by  a constant  c on  the  Interval  [s,S)  . To 
fix  the  value  of  c we  normalize  the  approximated  distribution 
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(11)  /^c  dx  + Pr{(S-y)  » S}  - 1 . 

Using  (10)  we  find  that 

Pr{(S-y)  - S}  - [1+M(D)]"^  , 
and  we  solve  (11)  for  c , yielding 

(12)  c - (1  - (1+M(D)]"^}/D  . 

Finally,  we  use  (9)  in  (12)  to  obtain 

(13)  c = m(y)/[l+M(D)]  ‘ (l-p)/D  , y e (o,D]  . 

We  use  (9)  and  (13)  in  (8)  to  get 

H = h{  l(l-p)/D]  ^(S-y-x)(t>  ^^^^^(x)dxdy 

+ p/^(S-x)(ti  ^^^^^(x)dx} 

(14)  B i p[H/h  + (L+l)p  - S + (l-p)D/2] 

pip  - [(l-p)/D]  /°  4'*^^'^^^(S-y)dy 

R i pK  . 

Finally  we  approximate  the  demand  distribution  with  a gamma  distribu- 
tion having  density  g(*  |u,6)  . Let 

x“  ^exp(-x/e)/ [1(0)6°*]  , x>0 

n . x<0 

^ (L+l)^^j  ^ G(x|a,6)  ‘ /!^„.g(y|Qi.B)dy  . 


4i*(L'’’1)  (x)  • g(x)a,6)  = 
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where 


a (L+l)u^/o  ^ 

B o2/y  . 

We  define  the  notation 

{f(x)  I*’  - f(b)  - f(a) 
a 

and  use  (15)  in  (14)  to  yield 

H = ph[SG(sla.B)  - aBG(s|(i+l,6)] 

+ [h(l-p)/2D]  {x^G(x|a,6)  - 2a BxG (x | a+1 , B) 
+ (a+l)aB^  G(x|a+2, B) 

(16) 

B i p[H/h  + (L+1)  - S + (l-f)D/2] 

P = Gisla.S)  - l(l-p)/D]{xG(x|u,B) 

- aBG(x|a+l,6)  1^ 

R « pK  . 


Observe  that  the  approximations  (L6 ) depend  on  the  policy  parameters, 
economic  parameters,  and  the  mean  and  variance  of  demand.  The  G-f unctions 
must  be  calculated  by  a numerical  procedure.  We  use  a series  expansion  for 
G(x|u,8)  when  x is  less  than  1 or  u6  , and  a continued-fraction 
expansion  otherwise.  The  procedure  Is  part  of  a package  of  computer 
programs  entitled  "The  IMSL  Library"  which  Is  marketed  by  the  International 
Mathematical  and  Statistical  Libraries,  Inc.,  Houston,  Texas. 


12 


Despite  the  effort  required  to  compute  G , the  expressions  In  (16) 
are  an  enormous  computational  simplification  over  (8).  In  Section  2.8 
we  discuss  the  possibility  of  using  a Normal  distribution  function  in  lieu 
ol  the  G-funcLion;  employing  the  Normal  distribution  would  facllliale 
manual  computations  of  the  approximations  we  derive  below. 

In  the  following  sections  we  use  expressions  (16)  to  develop 
regression  models  for  the  operating  characteristics.  We  fit  the 
parameters  of  the  regression  models  to  the  observed  characteristics 
of  576  items  (Table  2.1).  The  approximations  we  obtain  are  labeled 
with  subscript  "a"  when  they  are  used  for  the  entire  576-ltem  system. 
Subscripts  "a,p"  and  "a,o"  are  used  to  label  expressions  for  Power 
Approximation,  and  optimal  policies,  respectively. 

2 . 2 An  Approximation  for  ^e^^l e n_l_8hm^n_^  Cost 

We  use  (9)  in  (16)  to  obtain  the  expression  for  replenishment 

cost 

R i uK/lD  + (m4c2/u)/2]  . 

We  manipulate  the  expression  to  form  a linear  regression  model 

(mK/R)  - Aq  + A^  D + A^  P + A^  (a2/ii)  + E , 

where  are  constants  to  be  fit  and  l Is  the  error  term.  We 

use  least-squares  regression  to  fit  the  model  to  the  system  of  576 
inventory  policies  In  Table  2.1.  That  Is,  for  each  of  these  policies  we 
use  as  data  the  actual  values  of  pK/R  , D , p , and  o^/p  . The  result 
Is  the  following  numerical  approximation  for  R 
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(17)  R - Ku/(1.003D+.4942y+.499a)^/w-.5339)  , 
a 

which  has  a coefficient  of  determination  (fraction  of  variance  explained) 
of  0.9999  for  the  quantity  yK/R  . 


To  characterize  the  quality  of  the  approximation  we  define 
(18)  A(R  ) r 100%  X 1r  -r1/R  . 

In  other  words,  A(R  ) is  a percentage  measure  of  the  approximation 

a 

error  for  a particular  item.  The  system  of  576  items  yields  an  average 

value  for  A(R  ) of  0.1%.  Table  2.2  gives  the  distribution  of  values 
a 

for  A(R  ) . The  approximation  is  clearly  an  accurate  one  for  the  system, 
a 


Table  2.2 

Frequencies  of  A(R  ) in  a 576-item  System 


Range 

for  A(R  ) 
a 

Number 
of  Items 

Cumulative  Percentage 
of  Items 

[0%,2%) 

574 

100'% 

(2%, 4%) 

2 

100% 

The  largest  value  of  A(R^)  is  2.5%  and  it  is  attained  by  two  optimal 

policies:  each  has  a varlance-to  mean  ratio  of  1,  leadtiroe  0,  mean 
demand  16  and  setup  cost  ratio  32.  The  unit  penalty  cost  ratios  are 
24  and  99,  yielding  the  optimal  policies  (16,39)  and  (20,43),  respectively. 
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2 

Since  the  coefficients  of  D,  u,  and  o /p  in  (17)  are  so  close  to 
the  theoretical  values  of  1,  1/2,  and  1/2,  respectively,  we  consider  the 
following  regression  model. 

[pK/R  - D - l/2(p  + o^/p)]  - A + c, 

where  A is  a single  constant  to  be  fit  and  e is  the  error  term.  Using 
least-squares  regression  to  fit  the  model  to  our  576-ltem  system  we  obtain 
the  approximation 

r:  - Kp/[D  + (p  + o^/p)/2  - .5121]  . 

a 

Approximation  R"  provides  essentially  the  same  accuracy  as  (17).  The 

system  of  576  items  yields  an  average  value  of  A(R^)  of  0.1%  and  the  dls- 

a 

trlbutlon  of  values  for  A(R')  is  the  same  as  that  given  in  Table  2.2  for 

a ” 

A(R  ). 
a 


Ua 


2 . 3 An  Approximation  for  Holding  Cost 


We  can  treat  the  unit  holding  cost  as  a redundant  (normalizing) 
parameter  In  our  model,  and  so  we  divide  the  holding  cost  expression  In 
( 16)  by  h yielding 

H/h  = p[SG(S|u,ti)  - <iB  G(s|a+l,ti)] 

+ [(l-p)/2D]{x^G(x|a,ti)  - 2uBxG(x|a+l,  6) 

+ (a+l)u6^G(xl!i+2,6)  . 

We  take  advantage  of  our  improved  estimate  of  replenishment  frequency 
frf>m  (17)  and  replace  p with 

(19)  r R /K=p/(i.003D+.49A2u+.A99a' ^/u-.5339)  . 

a 

Tlie  result  is  a quantity  that  we  denote  as  W , given  by 

(20)  W r(SG(s|u,e)  - uBG(s|'i+l,e)J 

+ t(l-r)/2D]{x^G(x|a,B)  - 2aBxG(x|u+l. j) 

+ ('i+l)aP^  G(x|u+2,B)  Ig  . 

We  calculated  values  of  W in  the  576-ltem  system  and  in  comparing  them 
with  the  actual  values  of  H/h  , we  found  a systematic  variation  with 
respect  to  M and  0^/p  . This  motivates  the  linear  regression  model 

H/h  - Aq  + A^  W + A2  U + (o^/u)  + £ , 

where  Aq,...,A2  are  constants  to  be  fit  and  c Is  the  error  term. 

As  in  Section  2.2,  we  use  least-squares  regression  to  fit  the  model  to 
the  system  of  576  items.  The  result  Is  a coefficient  of  determination 
of  0.9999  for  the  approximation 
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0 


(21) 


I 


i 


H - h(W  -.1512  u + .1684o2/u  + .0689)  . 


a 


We  define  A(H  ) as  in  (18).  The  system  of  576  items  yields  an  average 
a 

value  of  0.7%  for  A(H  ) . The  distribution  of  values  for  A(H  ) is 

a a 

summarized  in  Table  2.3.  Our  approximation  appears  to  be  an  excellent 

one . 


Table  2.3 

Frequencies  of  A(H  ) in  a 576-item  System 




Range 

for  A(H  ) 
a 

Number 
of  Items 

Cumulative  Percentage 
of  Items 

[0%,2%) 

553 

96% 

[2%, 4%) 

17 

99% 

[4%, 6%) 

5 

100"% 

[6%, 8%) 

0 

ioo"% 

[8%, 10%) 

1 

100% 

The  largest  value  of  A(H  ) is  9.2%,  attained  by  an  optimal  policy  with 

a 

variance-to-mean  ratio  9,  leadtime  0,  mean  demand  2,  unit  penalty  cost 
ratio  4,  and  setup  cost  ratio  32.  In  geneial,  the  largest  errors  appear 
to  occur  for  high  values  of  variance-to-mean  ratio  and  low  values  of  all 
other  parameters. 


2.4  An  Approximation  for  Backlog  Protection 

Backlog  protection  is  the  frequency  of  periods  without  backlogs,  that 
is,  one  minus  the  backlog  frequency.  Since  it  is  a critical  measure  of 
service,  it  is  of  central  Interest  to  the  Inventory  systems  designer. 
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I III  Section  2.1  we  obtained  the  following  approximation  for  backlog 

^ protection 

P = (:(S|a,S)  - [(l-p)/D]  {SG(s]a,6)  - sG(s|a,B) 

- 'xBG(S  |a+l  ,6)  + (jBG  (s  I u+1 , 6)  } . 

We  replace  o with  r , as  given  by  (19),  to  obtain 

P = V . rG(s|a,B)  - [(l-r)/D]  (SG(s|u,B)  - sG(s|a,B) 

- 3BG(s|a+l,6)  + aiBG(s  |ai+l,B) } . 

We  form  a linear  regression  model  by  adding  terms  for  mean  demand  and  the 
var iance-to-mean  ratio  of  demand,  as  we  did  in  our  approximation  for 
I holding  cost 

(22)  P = Aq  + V + A2  p + (o2/p)  + E . 

This  model  yields  a coefficient  of  determination  of  only  0.669  when  fit 
to  the  576-item  system.  We  also  fit  (22)  separately  for  Power  Approxima- 
tion and  optimal  policies,  yielding  coefficients  of  variation  of  0.660 
and  0.681,  respectively.  These  poor  results  forced  us  to  seek  an 
alternative  approach. 

We  can  revise  the  regression  model  by  using  a theoretical  result. 

When  demand  is  continuously  distributed,  an  optimal  policy  yields 
(p/h)  / (1+p/h)  for  backlog  protection.  When  the  demand  distributions 
are  discrete,  (p/h)  / (1+p/h)  is  a lower  bound  on  P for  optimal  policies. 
It  was  observed  in  F.hrhardt  (1976)  that  the  Power  Approximation  and  op- 
timal policies  differed  in  tlieir  backlog  frequency  performance.  Therefore 
we  decided  to  fit  the  two  policy  rules  separately. 
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I 

We  develop  the  revised  linear  regression  model  by  multiplying  (22) 
by  (1+p/h)  to  yield 

(1+p/h)  P - Aq  + (p/h)  + A2  V + Aj  (Vp/h)  + A^  W 

+ Aj  (up/h)  + A^  (o2/u)  + A^  [(o^/u)(p/h))  + e , 

This  model  dramatically  Improves  the  fit  after  only  one  step  of  a 
stepwise-regression  computer  program.  For  optimal  policies,  the  simple 
expression 

(23)  '’a  o * (0-0857  + p/h)  / (1-Kp/h) 

yields  a coefficient  of  determination  of  0.99999  for  (l+p/h)P  . We 
have  the  same  coefficient  of  determination  for  Power  Approximation  policies 
wl  th 

(2A)  P « (0.0695  + p/h)  / (1-t-p/h)  . 

a , p 

We  define  ^(P^)  ^s  the  absolute  value  of  the  percentage  error  of 

approximation  when  (23)  and  (2A)  are  used.  Table  2.4  is  a summary  of  the 

distribution  of  ^(^g)  the  576-item  system.  The  average  value  of 

A(P  ) in  the  system  is  0.7Z.  All  items  with  A(P  ) larger  than  4X  have 
a a 

Power  Approximation  policies  with  p/h  equal  to  4.  The  approximation 
(23)  Is  somewhat  better  than  that  in  (24)  as  40  of  the  items  wltli  devia- 
tions larger  than  2%  are  for  the  Power  Approximation.  All  the  items 
with  A(P  ) larger  than  4%  have  Power  Approximation  policies  with  p/h 
equal  to  4.  Of  those  items  with  deviations  larger  than  2X,  7 have  p /li 
equal  to  9 and  41  have  p/h  equal  to  4.  Our  approximations  are  especially 
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accurate  for  large  p/h  . 


Table  2.4 

Frequencies  of  A(P  ) in  a 576-ltem  System 
£ 

Number 
of  Items 

528 
39 
8 
1 

The  errors  of  approximations  (23)  and  (24)  display  a mild  dependence 
upon  the  parameter  K/h.  We  therefore  seek  an  approximation  of  the  form 

P - (a^  + ajh/K  + p/h)/(l  + p/h)  , 

which  suggests  the  regression  model 

Y i [P  - p/(h  + p)]  = Aq  + A^(h/K)  + L. 

The  model  Is  fit  to  minimize  relat ive  squared  error  by  first  dividing  by 
the  left  hand  side,  yielding 

1 - Aq(1/Y)  + A^(h/KY)  + c/Y  . 

The  transformed  model  is  then  fit  using  least-squares  regression;  i.e., 
the  dependent  variable  Is  set  equal  to  one  for  each  item  and  the  Intercept 
Is  fixed  at  zero. 

We  obtain  the  following  approxlimat ions  for  optimal  and  Power  Approxi- 
mation policies,  respectively 


Range 

for  A(P  ) 
a 

[0Z.2X) 

(2Z.4Z) 

(4Z,6Z) 

(6Z.8Z) 
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0 " (-0543  + .9958  h/K  + p/h)/(l  + p/h) 

P'  - (-.0737  + 5.851  h/K  + p/h)/(l  + p/h)  . 
a , p 

We  define  the  absoulte  value  of  the  percentage  error  of  approxi- 

mation. The  average  value  of  A(P')  In  the  576-ltem  system  is  0.6X.  The 

di 

distribution  of  A(P^)  is  summarized  In  Table  2.4a.  Although  approximation 

3 

?'  does  not  appear  significantly  better  than  P on  the  average,  it  reduces 

3 3 

Table  2.4a 

Frequencies  of  A(P')  in  a 576-ltem  System 


Range 

for  (P  ) 
a 

Number  of 

1 1 ems 

Cumulative  Percentage 

of  Items 

[0%,  2%) 

529 

92% 

[2%,  4%) 

42 

99% 

[4%.  6%) 

5 

100% 

[6%,  8%) 

0 

100% 

the  number  of  items  with  errors  in  excess  of  4%  from  9 to  5.  We  recommend 
apr)ro:;i.;-.uitioii  P'  for  Power  Approximation  items  with  small  values  of  p/h(s9). 

2 5 An  Approximation  for  To t a 1 C^8_t 

We  obtain  an  expression  for  total  cost  by  summing  cost  components 
H,  B,  and  R,  and  using  approximations  (16)  for  B and  R 

T - H + B + R 

i (1  + p/h)  H + p((L+l)y  - S + (1-p)  D/2]  + pK  . 
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We  divide  by  h , replace  p with  r , as  given  by  (19),  and  use 
approximation  (20)  for  H to  obtain 

(25)  T/h  i (l+p/h)W  + p/h  ((L+l)p  - S + (l-r)D/2]  + rK/h  . 

As  we  discovered  in  obtaining  a fit  for  holding  cost,  a group  of  related 
terras  should  be  added  to  (25)  to  obtain  a good  fit  to  the  system's  data. 
The  linear  regression  model  we  employed  is 
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0 

T/h  » Aq  + W + A2  (Wp/h)  + Aj  ((L+l)yp/h]  + A^  (Sp/h) 
+ Aj  (Dp/h)  + A^  (rDp/h)  + A^  (rK/h)  + Ag  (p/h) 

+ A^  (rp/h)  + A^q  I(L+l)p]  + A^j  S + A^2  ® 

^ ^14  ^ ^ ^16  ■*■ 

+ Ajg[ (o^/p) (p/h) ] + e . 


We  fit  Che  model  to  the  system  of  576  items  using  stepwise  least- 
squares  regression.  The  following  expression  yields  a coefficient  of 
determination  of  0.998 

(26)  T - 1.110  hW  - .001049  pW  + .3364  Kr 

a 

- .2234  h + .3274  hD  + .4476  h o^/p  + .003062  p o^/p  . 

We  define  A(T  ) as  in  (18),  namely,  the  absolute  value  of  the  per- 
a 

centage  difference  between  and  T . The  average  value  of 

for  the  576-item  system  is  1.9%.  A summary  of  the  distribution  of 
A(T^)  for  the  system  is  given  in  Table  2.5.  The  fit  is  rather  good  with 

virtually  all  of  the  large  errors  occurring  for  large  values  of  a^/p  , 

and  small  values  of  L,  p,  and  p/h  . Of  the  63  items  having 

greater  than  4%,  52  are  for  the  Power  Approximation.  We  list  the  four 

items  with  A(T  ) exceeding  10%  in  Table  2.6. 
a 


Table  2.5 

Frequencies  of  A(T  ) In  a 576-ltem  System 


Range 

for  A(T  ) 
a 

Number 
of  Items 

Cumulative  Percentage 
of  Items 

[0Z.2Z) 

383 

66% 

(2Z,A2) 

130 

89% 

14Z.6%) 

45 

97% 

(6Z,K%) 

13 

99% 

(SZ.IOZ) 

1 

99% 

[10%.12Z) 

2 

100"% 

[12%. 14%) 

0 

100"% 

114Z,16Z) 

1 

ioo“% 

[162,18%) 

0 

ioo"% 

[18%. 20%) 

1 

100% 

Table  2.6 

Items  with  A(T  ) exceeding  lOZ 


Policy 

o^/b 

Power  Approximation 

9 

Power  Approximation 

9 

Power  Approximation 

9 

Power  Approximation 

1 

9 

K/h 

Error 

64 

10.3% 

64 

11.5% 

32 

14.6% 

32 

18.2% 

The  parameter  settings  contributing  most  often  to  large  errors  in 


(26)  are  Power  Approximation  Items  with  p = 2 and  p/h  » 4.  We  deleted 
all  18  items  having  their  combination  of  parameter  values  and  refit  (26) 
to  obtain 

T'  = 1.111  hw  - .00105  pw  + .3370  Kr 
a 

- .2190  h + .3269  hD  + .4775  hu^/p  + .00306  po^/p  , 

which  has  an  average  value  of  ^(T^)  of  1.8%  in  the  558-ltem  system.  A 
suranuiry  of  tlie  distribution  of  A(T')  is  given  in  Table  2.6a. 

Table  2.6a 

Frequencies  of  A(T^)  in  a 558-item  System 

3 


Range 

for  A(t') 
a 

Number 
of  Items 

Cumulative  Percentage 
of  Items 

[0%,  2%) 

377 

68% 

[2%,  4%) 

126 

90% 

[4%,  6%) 

42 

98% 

[6%,  8%) 

11 

100% 

[8%,  10%) 

0 

100% 

[10%,  12%) 

1 

100% 

[12%,  14%) 

0 

100% 

[14%,  16%) 

1 

100% 

Although  the  approximation  appears  to  be  accurate,  we  believe  that 
it  will  degrade  for  significantly  non-optlmal  policies.  The  differences 
between  (25)  and  (26)  and  the  pattern  of  values  for  A(T  ) suggest  that 

the  economics  of  optimal  policies  is  intrinsic  to  the  approximation 
obtained.  The  robustness  of  (26)  is  discussed  explicitly  in  Section  2.7. 

2.6  An  Approximation  Xpr  Backlog  Cjost 

In  Section  2.1  we  obtained  the  approximation  for  backlog  cost 

B - p[H/h  + (L+l)y  - S + (l-p)D/2]  . 

Using  the  approach  of  Section  2.3,  we  obtain  the  linear  regression  model 
B/h  • (Wp/h)  + A^  [ (L+Dyp/hJ  + A^  (Sp/h) 

+ A^  (yp/h)  + A^  [ (o^/y ) (p/h) ] + A^  (Dp/h) 

+ A^  (rDp/h)  + t , 

where  r is  given  by  (19),  and  e is  the  error  term.  When  fit  to  the 
system  of  576  items  the  result  is  a low  coefficient  of  determination  of 
0,444.  Splitting  the  576-ltem  system  into  a Power  Approximation  system 
and  an  optimal  policy  system,  and  then  fitting  separately  for  each  value 
of  p/h  , does  not  yield  a significant  Improvement.  The  two  systems 
also  were  fit  separately  for  each  value  of  o^/y  without  success. 

We  obtain  a better  fit  by  employing  the  identity 
B - T - H - R . 

Our  first  attempt  used  (17),  (21),  and  (26)  in 

B'  - T - U - R . 
a a a 
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Unf  orlunati'ly , B'  deviates  from  B by  an  average  of  18%. 

We  form,  instead,  a linear  regression  model  for  B using  nil  the 
variables  appearing  In  our  models  for  T,  H,  and  R 

(27)  «/b  “ Ay  + W + A^  (Wp/h)  + A^  [(L+l)pp/h]  + A^  (Sp/h) 

+ A^  (I)p/h)  + A^  (rDp/h)  + A^  (rK/li)  + A^  (p/h) 

+ Ag  (rp/h)  + Ajy  [(L+l)pJ  + A^^  S + A^^ 

A^^  (Dr)  + Aj^^  r + A^^^  p + A^^  (o^/p) 

+ A^^  (pp/h)  + A^g  I (p/h) (o^/p)]  + c . 


When  fit  to  the  576-item  system,  (27)  yields  a coefficient  of  determination 
of  0.957.  To  obtain  a reasonably  accurate  approximation,  we  further  split 
the  data  set  Into  Power  Approximation  and  optimal  policies,  and  fit  (27) 
separately  for  each  value  of  p/h  . 

We  obtain  average  coefficients  of  determination  of  0.999,  and  0.998 
for  optimal,  and  Power  Approximation  policies,  respectively.  The  fits 
are  of  the  form 

(28)  B /h>»  a,  W + a„  (L+l)p  + a^  S + a.  D + ar  Dr 

a 1 2 J 4 J 

+ a^  r + a^  Kr  + ag  p + a^  o^/p  + . 

We  list  the  fitted  values  of  {a^,l"l,10}  in  Table  2.7a  for  each  of  the 
eight  regressions. 
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Table  2.7a 


Coeff  icient  Values  for  Backlog  Cost  Expression  (2ti) 


Policy 

p/ll 

'*1 

^2 

■"3 

“5 

"7 

*8 

‘‘9 

^0 

Optimal 

3.231 

3.212 

-3.210 

1.680 

-1.809 

.1968 

-.1658 

.2600 

-.0127 

-.3447 

7.452 

7.473 

-7.470 

3.854 

-4.046 

.4851 

.2325 

. 3638 

-.0437 

-.4972 

24 

17.39 

17.37 

-17.37 

8.788 

-9.026 

1.034 

-.2005 

.3445 

.0336 

-.7617 

99 

62.86 

62.85 

-62.84 

31.56 

-31.74 

.9851 

-.2572 

.3412 

.0927 

-1.084 

Power 

4 

4.248 

4.276 

-4.278 

2.141 

-1.740 

-1.272 

-.0758 

-.1326 

-.1099 

-.0900 

Approx . 

9 

8.214 

8.192 

-8.194 

4.152 

-3.813 

-.3054 

-.1565 

-.1044 

-.0635 

-.4659 

24 

19.18 

19.14 

-19.14 

9.635 

-9.228 

-.3749 

-.1844 

-.1606 

-.0059 

-.6567 

99 

70.16 

70.18 

35.22 

-34.55 

1.428 

-.2483 

-.3159 

.1841 

-1.327 

As  the  high  coefficients  of  determination  indicate,  the  fits  are 

excellent  in  terms  of  absolute  errors  of  approximation.  To  Illustrate, 

the  average  absolute  error  of  approximation  for  288  optimal  policies  Is 

0.17  as  compared  with  the  corresponding  average  backlog  cost  5.69. 

Nevertheless,  because  many  of  the  backlog  costs  are  numerically  small 

values  and  the  magnitude  of  the  approximation  errors  does  not  seem  to 

depend  strongly  on  the  magnitude  of  the  backlog  cost,  the  approximation 

displays  some  ’ ^rge  relative  errors.  We  define  (18) » the 

absolute  value  of  the  percentage  difference  between  and  B for  a 

given  Item,  Tlie  average  value  of  A(B  ) in  the  576-ltem  system  is  4.  IX; 

more  specifically,  the  average  is  4.3%  for  optimal  policies  and  3.9%  for 

Power  Approximation  policies.  We  summarize  the  distribution  of  A(B  ) 

a 

2A 


a 


in  Table  2.8.  The  largest  values  of  occur  when  o^/y  equals  1 

and  p/h  equals  99.  All  16  Items  having  errors  in  excess  of  207.  have 
o^/u  equal  to  1,  and  13  of  them  have  p/h  equal  to  99.  When  we  consider 
only  those  items  with  o^/y  greater  than  1,  the  average  value  of  A(B  ) 
drops  to  2.5%,  but  there  are  still  11  out  of  384  items  with  deviations 
exceeding  10%.  When  we  consider  only  those  items  with  p/h  less  than  99, 
tile  average  value  of  drops  to  2.4%. 

To  improve  on  the  approximation  when  p/h  equals  99,  we  exhibit  in 

Table  2.7b  an  approxlntition  for  this  value  of  the  unit  penalty  cost  ratio 

that  also  depends  on  0^/y  . The  average  value  of  A(B  ) for  this 

a 

approximation  is  4.2%,  with  16  of  the  144  items  having  deviations  in 
excess  of  10%. 

We  attempted  to  ' . r .i.iation  by  fitting  to  minimize 

relative  squared  error  is  in  Section  2.4  (p.l9).  The  resulting  approxima- 
tion was  less  accurate  than  those  reported  above. 

Since  this  approxinva tlon  is  somewhat  less  accurate  than  the  others 
we  have  derived,  we  note,  that  we  liave  a good  alternate  measure  of  backlog 
performance  in  our  approximation  for  P , (23)  and  (24). 


Table  2.7b 

Coefficient  Values  tor  Backlog  Cost  Kxpression  (28) 
with  p/h  ° 99 


Policy 

o^/y 

^2 

— 

■'*3 

B 

— 

D 

^8 

“9 

^10 

Optimal 

1 

71.44 

71.37 

-71.38 

35.76 

-35.99 

2.246 

-.1238 

.2124 

-.4770 

-1.084 

3 

83.12 

83.15 

-83.14 

41.63 

-41.96 

3.309 

-.0362 

.2961 

-.1453 

-1.084 

9 

75.13 

75.14 

-75.13 

37.62 

-37.50 

1.764 

-.1130 

.0516 

.0657 

-1.084 

Power 

1 

59.44 

59.42 

-59.41 

29.71 

-29.23 

1.090 

-.0428 

-.3647 

.4862 

-1.327 

Approx. 

3 

80.31 

80.34 

-80.33 

40.21 

-4U.02 

2.361 

-.0853 

-.0832 

.0467 

-1.327 

9 

86.35 

86.34 

-86.35 

43.20 

-42.69 

1 1.326 

-.0577 

-.2866 

-.0224 

-1.327  j 

r 

i 


Number 
of  Items 


Cumulative  Percentage 
of  Items 


Range  tor 

A(B  ) 
a 


[0%,22:) 

275 

48% 

127:.  4%) 

147 

73% 

[4%,  67!) 

51 

82% 

[6%, 8%) 

27 

87% 

[8%. 10%) 

18 

90% 

[10%, 20%) 

42 

97% 

[20%, 30%) 

8 

99% 

[30%, 40%) 

6 

ioo"% 

[40%, 50%) 

1 

ioo’% 

[50%, 60%) 

0 

iOO‘% 

[60%, 70%) 

0 

ioo"% 

[70%, 80%) 

1 

100% 

2 . 7 Measures  of  Robustness 

We  summarize  our  results  by  listing  the  accurate  approximations  derived 
in  this  section. 

H - hW  - .1  2 hp  + .1684  ha^/p  + .0689  h 

a 

P » (0.0857  + p/h)/(H-p/h) 

a ,o 

P - (0.0695  + p/h)/(l+p/h) 

a , p 

R^  - Kp/(1.003D  + .4942P  + .4990  o^/p  - .5339) 

B • h[a,  W + a-  (L+l)y  + a,.  S + a.  D + a.  Dr  + r 

a 12  J 4 3 0 

+ a^  Kr  + ay  p + O^/p  + 
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t 


where 


hd.llO  W + .327A  D + .4476  o^/u  - .2234) 

+ p(-.001049W  + .003062  a^/p)  + .3364  rK  , 


k * 


T 


a 


► 

I 


r = R / K 
a 

W = r[sr,(s|a,B)  - aBG (S  | a+1 , 6)  ] 

+ [U-r)/2D]  {x^G(x|fi,S)  - 2oiBxG  (x  | a+l , g) 
+ (a+l)ni6^G(x  |a+2,  B) 


G(x|(j,B) 


I [r(a)B'*]  ^ ^exp(-y/6)dy 


\ 0 

a = (L+Dp^/o^ 

B = a^/u  , 


X > o 


X < o 


and  the  coefficients  for  B are  given  in  Tables  2.7a  and  2.7b. 

a 

Recall  that  H , R , and  T can  be  •’"piled  to  both  optimal  and  Power 
a a a 

Approximation  policies. 

In  Tables  2.9  and  2.10  we  summarize  the  multi-item  system  accuracy 
of  our  approximations  as  a function  of  the  parameter  settings  (see 
Table  2.1).  Each  entry  Is  the  percentage  excess  of  the  sum  of  approx- 
imate values  over  the  sum  of  actual  values  for  all  the  items  having  the 
given  parameter  setting.  The  approximation  we  use  for  backlog  cost 
includes  the  Improved  fit  for  p/h  e(|ual  99  (see  Table  2.7b).  Observ*' 
that  when  the  approximations  are  aggregatc’d  over  several  parameter  settings 
that  the  errors  cancel  to  a large  extent  so  that  most  of  the  relative 
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Percentage  Errors  of  Analytic  Approximations  for  a 288-iten  System 


errors  in  these  tables  are  well  under  1%.  Even  for  the  backlog  cost 


approximation,  the  relative  errors  for  multi-item  aggregates  are 
tolerable. 

We  test  the  robustness  of  the  approximations  by  using  them  in  a 
multi-item  system  with  the  parameter  settings  of  Table  2.11.  Note  that 
all  the  numerical  parameters  have  values  not  found  in  the  576-ltem  system. 


Table  2.11 

A b4-item  System  wl Lh_  New  Parameter  Settings 


Factor 

Levels 

Number  of  Levels 

Dem.and  distribution 

Negative  Binomial  (o^/y  •=  5) 
Negative  Binomial  (o^/y  = 15) 

2 

Mean  demand 

0.5,  7.0 

2 

Replenishment  lead time 

1,  6 

2 

Replenishment  setup  cost 

16,  48 

2 

Unit  penalty  cost 

49,  132 

2 

Unit  holding  cost 

1 

1 

Policy 

Optimal  policy. 

Power  Approximation  policy 

2 

Each  parameter  has  one  Interpolated  value  and  one  extrapolated  value. 

A full  factorial  combination  of  the  values  is  used,  yielding  64  items. 
The  system  is  a rather  severe  test  of  robustness  since  only  two  items 
Ivave  all  parameters  with  values  within  the  ranges  used  to  derive  the 
approximations.  There  are  10  items  with  one  extrapolated  parameter, 

20  items  with  two  extrapolated  parameters,  20  with  three  extrapolations, 
10  with  four  extrapolations,  and  2 items  with  all  five  parameters  extra- 
polated . 

iO 


We  compare  actual  values  of  H,  P,  R,  and  T for  the  64  Items  with 
our  analytic  approximations.  Backlog  cost  B Is  not  considered  because 
of  the  complexity  of  our  approximation  and  the  absence  of  an  explicit 
dependence  on  p/h  In  (28).  The  average  percent  deviations  from  actual 
values  of  H,  P,  R,  and  T are  1.6%,  0.2%,  1,4%,  and  2,6%,  respectively. 

The  distributions  of  percent  deviations  are  summarized  In  Table  2,12.  Our 
approximations  are  quite  accurate  considering  the  wide  range  of  parameters 
spanned  by  the  system. 

The  holding  cost  approximation  Is  extremely  accurate  for  all  cases 

with  u greater  than  0.5  or  o^/p  less  than  15.  All  items  with  deviations 

greater  than  4%  liave  ; equal  0.5  and  o^/p  equal  15.  If  we  consider 

only  the  items  with  fewer  than  two  parameters  extrapolated,  the  average 

A(H  ) Is  0.4%. 
a 

Table  2.12 

Percentage  Deviations  of  Analytic  Approximations 
in  a 64-ltem  System 


Range  of 
Deviation 

Holding 

Cost 

Backlog 

Protection 

Replenishment 

Cost 

Total  Cost 

(0%,2%) 

48  (75%) 

64  (100%) 

48  (75%) 

30  (47%) 

12%, 4%) 

6 (84%) 

8 (88%) 

22  (81%) 

[4%, 6%) 

5 (92%) 

o 

00 

00 

6 (91%) 

1 6%, 8%) 

3 (97%) 

6 (97%) 

4 (97%) 

(8%, 10%) 

2 (100%) 

2 (100%) 

1 (98%) 

[10%, 12%) 

1 (100%) 
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I 

The  backlog  protection  approximation  is  excellent,  with  only  one 
item  having  a deviation  in  excess  of  0.7%. 

Our  approximation  for  replenishment  cost  is  also  robust.  All  items 
with  deviations  In  excess  of  4%  have  p equal  0.5,  0^/p  equal  15,  and 
K/h  equal  lb.  Items  with  fewer  than  two  extrapolated  parameters  have 
an  average  of  0.1%. 

Low  p and  high  o^/p  are  also  sources  of  large  errors  for  our 
total  cost  approximation.  All  items  with  deviations  in  excess  of  4%  have 
I either  p equal  0.5  or  o^/p  equal  15,  or  both.  Items  with  fewer  than 

two  extrapolated  parameters  have  an  average  deviation  of  1.2%. 

We  commented  in  Section  2.5  that  the  approximation  for  total  cost 
i may  degrade  for  significantly  non-optlmal  policies.  The  remark,  is 

equally  valid  for  the  backlog  protection  expressions  (23)  and  (24),  since 
they  are  based  on  a theoretical  result  for  optimal  policies.  We  now 
proceed  to  examine  the  accuracy  of  the  approximations  for  non-optimai 
policies . 

Consider  the  following  system  of  non-optimal  policies.  We  choose 
a base-case  Item  with  o^/p  equal  5,  p equal  9,  L equal  2,  p/h  equal 
49,  and  K/h  equal  48.  The  optimal  policy  for  this  item  is  (43,73).  We 

now  use  this  policy  on  items  with  different  parameter  values.  The  new 
parameters  are  obtained  by  increasing  and  decreasing  each  basii-case 
parameter  value,  one  at  a time,  yielding  10  items;  the  para.Tie*.er  values 
of  the  system  are  displayed  in  Table  2.13.  For  each  Item  we  compare  the 
actual  and  approximate  values  of  H,  P,  R,  and  T. 
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Percentage  Errors  of  Approximations 


Changed 

Value 

0^/u  = 4 

(-20%) 

6 

(+20%) 

M = 7 

(-22%) 

11 

(+22%) 

L =1 

(-50%) 

3 

(+50%) 

p/h  = 39 

(-20%) 

59 

(+20%) 

K/h  "=  38 

(-21%) 

58 

(+21%) 

Average  of 

Absolute  Values 

Holding 

Cost 


Backlog 

Protection 

Replenishment 

Cost 

Total 

Cost 

-.6% 

-.00% 

6.0% 

.7% 

.00% 

-5.0% 

-1.3% 

-.03% 

13.7% 

2.9% 

.03% 

-22.2% 

-1.6% 

.00% 

12.6% 

5.5% 

.00% 

-36.2% 

-.5% 

.00% 

3.9% 

.3% 

.00% 

-2.2% 

.0% 

.00% 

4.0% 

.0% 

.00% 

-2.2% 

1.3% 

.01% 

- -1 

10.8% 

The  approximations  for  holding  cost  and  replenishment  cost  are  very 
accurate,  with  average  percentage  deviations  of  0.04X  and  0.01%, 


respectively.  The  approximation  for  backlog  protection  is  somewhat  leas 
accurate,  with  the  largest  errors  occurring  for  large  values  of  leadtime 
and  mean  demand.  The  total  cost  appioxlmation  does  not  perform  well  in 
the  system,  deviating  by  an  average  of  10.8%.  Thus  we  conclude  that  the 
approximations  for  backlog  protection  and  total  cost  should  be  used  with 
caution  for  significantly  non-optimal  policies.  An  approach  to  reducing 
the  errors  might  be  gleaned  from  the  pattern  of  deviations  in  Table  2.13. 
Notice  that  when  each  parameter  is  larger  than  in  the  base  case,  the 
approximation  underestimates  the  total  cost,  and  when  the  parameter  is 
smaller  than  in  the  base  case,  the  approximation  overestimates  the  total 
cost.  The  reverse  is  true  for  backlog  protection. 

2 . 8 Topi c ^_f or  Future  Research 

In  this  final  portion  of  Section  2 we  suggest  several  extensions  of 
our  research  on  analytic  approximations. 

Notice  in  the  summary  of  the  approximations  in  Section  2.7  that  the 
G-functlon  is  used  directly  only  in  the  calculation  of  W . As  we  mentioned 
in  Section  2.1,  the  computation  of  G is  complex,  even  on  a digital 
computer.  As  an  alternative  to  W for  manual  computations,  we  could 
calculate 

W = r[SN(S  ) - aeN(S,)] 
o 1 

+ ((l-r)/2D]  [S^N(S  ) - 2aBSN(S,)  + (a+l)aB^N(S„) 
o J ^ 

- 8^N(s  ) + 2«BsN(s,)  - (a+l)aB^N(s  )] 
o 1 2 

where  N(‘)  is  the  unit  Normal  distribution  function,  and 
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= [S  - (a+i)  ] / [(u+l)B^] 

s.  = [s  - („+i)  ] / [(a+l)6^]'^  . 

Fvirthe.r  research  is  needed  to  Investigate  how  much  degradation  in  the 
approximations  occurs  if  W'  is  used  instead  of  W . Also,  it  would 
be  possible  to  refit  the  coefficients  in  the  approximations  for  H,  B, 
and  T utilizing  W'  instead  of  W . 

We  now  discuss  an  approach  for  improving  the  approximations  for 
P and  B . Recall  that  these  fits  have  some  items  with  large  relative 

errors  despite  very  high  coefficients  of  determination.  That  is,  the 

regressions  were  very  effective  in  minimizing  absolute  squared  errors, 
but  several  items  with  small  values  of  P or  B were  left  with  large 
relative  errors.  Additional  research  is  needed  to  devise  a method  to  fit 
these  approximations  to  data  such  that  relative  squared  error  is  minimized 
A final  research  topic  involves  the  potential  use  of  the  approx- 
imations in  forecasting  the  future  performance  of  Inventory  systems.  The 
method  of  forecasting  discussed  in  F.hrhardt  (1976)  is  based  on  a retrospec 
tive  simulation  of  system  performance  using  a sequence  of  past  demands. 

In  actual  practice,  this  method  can  be  very  expensive  since  a time  series 
of  demand  data  must  be  stored  for  every  item  Involved  in  the  forecasting. 
If  the  approximations  of  this  chapter  are  substituted  for  the  retrospectlv 
simulation,  only  the  sample  mean  and  sample  variance  of  demand  must  be 
stored.  Additional  research  Is  needed  to  evaluate  the  accuracy  of  this 
method  of  forecasting. 
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i.  MDLTIPLICATIVK  APPROXIMATIONS  FOR  FIXED  AND  STATISTICAL  POLICIES 


In  Secf.Ioii  2 we  derived  analytically-motivated  approximations  for 
operating  characteristics  of  fixed  (s,S)  policies.  The  approximations 
are  accurate  but  they  lack  the  ease  of  interpretation  and  calculation 
that  we  desire  for  sensitivity  analysis  of  systems  design  Issues. 
Furthermore,  they  inherently  do  not  treat  statistically-derived  policies, 
because  they  are  based  on  theory  that  is  valid  only  foi  fixed  [)olicles. 

In  this  chapter  we  develop  approximate  expressions  for  the  operating 
characteristics  (7)  of  both  fixed  and  statistical  policies.  We  postulate 
a functional  form  for  the  expressions  from  examining  empirical  observa- 
tions and  derive  explicit  approximations  using  least-squares  regression. 
The  expressions  are  particularly  attractive  for  sensitivity  analysis 
because  they  are  products  of  simple-to-compute  functions,  each  of  which 
depends  on,  at  most,  two  parameters.  Further,  they  do  not  require  the 
calculation  of  the  replenishment  policy  itself. 

As  In  Section  2,  we  obtain  good  approximations  for  holding  cost, 
backlog  protection,  replenishment  cost,  and  total  cost.  On  the  average, 
the  fixed  policy  approximations  deviate  from  the  actual  values  of  these 
characteristics  by  4%,  1%,  2%,  and  3%,  respectively.  The  statistical 
policy  expressions  have  average  errors  of  4%,  1%,  2%,  and  3X,  respec- 
tively. Using  the  approach  in  this  chapter,  we  have  not  found  a good 
approximation  for  backlog  cost. 
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3 . 1 Methodology 

Our  motivation  here  derives  from  the  empirical  observations  found 
in  MacCormick  (1974),  Estey  and  Kaufman  (1975),  and  Ehrhardt  (1976). 

In  all  these  Instances,  multi-item  systems  were  studied  with  a 
full  factorial  setting  of  parameters,  and  constant  demand  variance-to- 
mean  ratio.  For  each  of  the  systems,  the  average  total  cost  of  an  item 
appears  approximately  to  follow  a multipllcatively  separable  function  of 
the  economic  and  demand  parameters.  The  relationship  has  been  observed 
when  the  systems  are  controlled  with  optimal  policies,  approximately 
optimal  policies,  and  statistically-derived  policies. 

For  an  item  controlled  with  an  optimal  or  approximately  optimal 
policy,  we  postulate  that  the  average  total  cost  per  period  can  be 
approximated  by  a function  of  the  form 

(29)  T/h  = cf ^ (L,02/p)f 2(p.o^/M)f2(p/h,a^/M)f^(K/h,o^/p)  . 

Expression  (29)  Is  consistent  with  the  empirically-observed  property  of 
multiplicative  separ.abl lity  for  a given  value  of  o^/p  . We  allow  o^/u 
to  enter  all  four  functions  because  we  have  no  reason  to  expect  separa- 
bility in  this  parameter.  For  subsequent  analysis,  we  limit  consideration 
to  functions  of  the  form 

Y (x,o^7p) 

(30)  - X exp[6^(x,o2/lj)] , 1 » 1,...,4  , 

where  y^(x,o^/p)  and  6^(x,o^/p)  are  linear  combinations  of  the  variables 

I,  X,  1/x,  o^/u,  p/a^  , 

(11)  x^,  1/x^,  (o^/p)^,  (p/o^)^  , 

xo^/p,  xp/a^,  o^/px,  p/o^x  . 
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As  In  Section  2,  we  transform  the  characteristics  (7)  to  H/h, 

B/h,  P(l+p/h),  R/h,  and  T/h  . We  seek  approximations  for  all  five 
(juantlties  using  linear  least-squares  regression  on  models  of  the  form 
(29),  (30),  (31).  Given  the  multiplicative  form  in  (29),  we  make  a 
logarithmic  transformation  before  applying  least-squarea  computations; 
hence,  in  reporting  below  the  resulting  coefficients  of  determination, 
we  always  shall  be  referring  to  the  proportion  of  variance  explained  for 
the  logarithmically  transformed  variables.  The  576-item  system  of 
Table  2.1  provides  the  data  source  for  separate  fits  for  the  character- 
istics of  optimal  and  Power  Approximation  policies. 

For  the  statistical  policy  characteristics  (see  Section  1 ) , we  use 
the  72-ltem  systems  of  Table  3.1.  We  simulated  the  system  with  a^/p 
equal  3 using  a 26-period  revision  Interval.  For  the  system  witli  o^/p 
equal  9,  we  utilized  revision  intervals  of  13,  26,  and  52  periods.  We 


Table  3.1 


72-item  System  Parameters 


Factor 

Levels 

Number  of  levels 
in  each  system 

Demand  distribution 

Negative  Binomial  (a^/p“9) 
or  Negative  Binomial  (a^/p“3) 

1 

Mean  demand 

2,  4,  8,  16 

4 

Unit  holding  cost 

1 

1 

Unit  backlog  penalty  cost 

4,  9,  99 

3 

Replenishment  setup  cost 

32.  64 

2 

Replenishment  leadtime 


0.  2,  4 


seek  approximations  for  statistical  policy  operating  characteristics  by 
performing  separate  fits  for  each  of  the  72-item  systems.  We  use  regres- 
sion models  of  the  form  (29),  ( JO)  with  and  (5^  being  linear 

combinations  of  the  variables 

* 

(32)  1,  X,  1/x,  x^,  1/x^  . 

For  several  operating  characteristics  we  seek  an  explicit  functional 
relationship  with  revision-interval  length  T . These  operating 
characteristics  are  fit  in  the  o^/p  equal  9 system  using  y^  and  6^ 

that  are  linear  combinations  of  all  the  variables 

(33)  1,  X,  1/x,  T,  1/T, 

2 2 
X , 1/x 

As  we  shall  observe,  several  fits  for  individual  y^  and  6^  do  not 
do  not  require  the  presence  of  all  these  variables. 

3 . 2 Approximations  for  Fixed  Policies 

In  this  section  we  describe  approximations  for  fixed  policy 
operating  characteristics.  Since  optimal  and  Power  Approximation 
policies  are  fit  separately,  we  adopt  the  convention  that  a subscript 
"o"  labels  an  optimal  policy  approximation,  and  a subscript  "p"  labels 
a Power  Approximation  quantity.  In  some  instances  we  use  summaries  that 
combine  results  for  optimal  and  Power  Approximation  systems.  We  use  the 
subscript  "f"  to  label  combined  fixed-policy  results. 

3.2.1  Average  Holding  Cost  £er_P^rio^ 

The  following  approximation  for  optimal  holding  cost  yields  a 
coefficient  of  determination  of  0.992  for  the  288  data  points  < 
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= 2.022  h exp(.09835o2/u)n’^^®®(L+l) 

(p/h)-2.A24../p-.04409y/a2+.01552o2/y 

^l^^^j.3660-.0269:  •■  /u 

On  the  average,  is  within  4.2%  of  H for  the  system. 

We  obtain  the  following  approximation  for  Power  Approximation 
holding  cost 

u o cQo  u / no/  2/  \ . 5076-.004964o2/p 
= 2.598  h exp(.1204o"^/u)p 

(L^l).2287+.01472o2/p-.0900p/a2 

^p^j^^-3.040h/p+.01137a2/M-.07292p/a2+.03355(a2h/pp) 
.3057-.02667o2/p+.03058u/o2 

» 

which  yields  a coefficient  of  determination  of  0.994  and  is  within  3.9% 
of  H , on  the  average. 

The  distribution  of  errors  for  the  optimal  and  Power  Approximation 
fits  are  so  alike  that  we  summarize  them  as  one  system.  Define  A(Hj) 

as  in  (18),  the  absolute  percentage  difference  between  or  and  H 

for  an  item.  Table  3.2  is  a summary  of  the  distribution  of  A(Hj)  in  the 
system.  The  largest  values  of  A(Hj.)  occur  consistently  for  a varlance- 
to-mean  ratio  9,  leadtime  of  0,  unit  setup  cost  ratio  of  32,  and  penalty 
cost  ratios  of  4 and  9.  This  is  a pattern  we  noted  also  in  Section  2.3 
for  the  analytic  approximation  of  holding  cost. 

We  attempted  to  Improve  the  approximation  by  refitting  separately  for 
each  value  of  leadtlmc  I..  The  resulting  approximations  were  not  signlf- 


(K/h) 


leant ly  better  than  IK  and  H . 

0 p 
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Table  3.2 

Frequencies  of  A(H^)  in  a 576-ltem  System 


Cumulative  Percentage 
of  Items 


28% 

60% 

77% 

87% 

94% 

98% 

99% 

99% 

100"% 

100"% 

100"% 

100~% 

100"% 

100"% 

100% 


3.2.2  Averaee  BackioK  Cost  per  Period 


We  have  not  found  a very  good  multiplicative  fit  for  backlog  rose. 
The  result  is  not  surprising  since  a multiplicative  form  Is  not  observed 
empirically  for  this  character Istlc . 

In  fitting  a multiplicative  model  of  the  form  (29)*  (30),  (31)  to 


backlog  cost,  we  obtained  coefficients  of  determination  of  0.976  and 

0.952,  respectively,  for  optimal  and  Power  Approximation  policies.  An 

improvement  resulted  when  we  ran  the  regressions  separately  for  each 

value  of  p/h.  We  obtained  average  coefficients  of  determination  of 

0.988  and  0.981  for  optimal  and  Power  Approximation  policies  with  a 

t unction  of  the  following  form 

b -b  o^/p  c.+c.o^/p+c.p/o^ 
a^/h  - a(a2/p)  ^2  ^12  3 

d +d  u/a2  e 02/p+e  p/o2 

(L+1)  (K/h) 

Values  of  the  coefficients  a,  bj^,...,e2  are  listed  in  Table  3.3  for 
each  policy  and  each  value  of  p/h  . 

Table  3.3 

Values  of  Coefficients  in 


Policy 

p/h 

B 

B 

^2 

"l 

^^2 

"3 

^1 

^^2 

r 

^1 

^2 

Optimal 

4 

. 3398 

1.003 

-.0100 

.5240 

-.0119 

.0056 

.4239 

-.1792 

-.0008 

.2891 

9 

.5800 

.7796 

.0034 

.4909 

-.0189 

-.0321 

.2999 

.0528 

.0034 

.1063 

24 

.2297 

1.706 

-.0350 

.3310 

-.0104 

.1854 

.2318 

.1799 

-.0033 

.2646 

99 

.3222 

-.0303 

.2786 

-.0104 

.1669 

.1739 

.2939 

-.0020 

.1378 

Power 

4 

.3704 

1.167 

-.0806 

.3775 

.0070 

-.0296 

.3174 

-.0074 

.0280 

.2748 

Approx . 

9 

.3316 

1.226 

-.0596 

.3711 

.0043 

.0168 

.3870 

-.0018 

.0162 

.2557 

24 

.2539 

1.520 

-.0456 

.3759 

-.0111 

.1135 

.3237 

.1619 

.0077 

.2296 

99 

. 1532 

2.235 

-.0875 

.2144 

-.0258 

.4537 

-.0614 

.7401 

.0240 

.1297 
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On  the  average,  approximation  deviates  from  actual  values  of 


B by  6.2%.  In  Table  3.4  we  summarize  the  distribution  of  percentage 
deviations  in  the  multi-item  system.  We  see  that  the  approximation  is 
not  accurate,  with  18%  of  the  deviations  exceeding  10%,  and  5%  of  the 
deviations  larger  than  20%.  The  pattern  of  errors  is  similar  to  that 


Table  3.4 

Frequencies  of  A(Bj)  In  a 576-ltem  System 


Range  for 

A(Bj) 

Number 
of  Items 

Cumulative  Percentage 
of  Items 

[0%,2%) 

164 

28% 

f2%,4%) 

116 

49% 

[4%, 6%) 

85 

63% 

[6%, 8%) 

53 

73% 

[8%, 10%) 

52 

82% 

[10%, 20%) 

78 

95% 

[20%, 30%) 

22 

99% 

[30%, 40%) 

5 

100'% 

[40%, 50%) 

0 

ioo"% 

[50%, 60%) 

1 

100% 

noted  In  Section  2.6  for  the  analytic  approximation  for  backlog  cost. 
The  largest  errors  occur  for  items  having  high  p/h  and  low  o^/u  . 
We  tried  improving  the  approximation  by  fitting  separately  for  each 
combination  of  values  of  p/h  and  a^/p  , but  the  results  did  not 
improve  significantly. 
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We  Improve  the  approximation  by  fitting  to  minimize  relative  squared 


error  as  in  Section  2.4  (p.l9).  We  obtain  a separate  fit  of  the  following 

2 

form  for  each  combination  of  values  of  p/h  and  o /p 

Bj/h  = ap'’(L+l)‘^  (K/h)‘'  . 

2 

The  result.s  did  not  Improve  for  Items  with  a =1  but  the  average  per- 

2 

centage  errors  dropped  to  5.7%  and  3.3%,  for  items  with  a /p  equal  to  3 

2 

and  9,  respectively.  Values  of  the  coefficients  a,b,c  and  d for  o /p 
equal  to  3 and  9 are  listed  in  Table  3.4a  and  the  distribution  of  percent- 
age errors  is  summarized  in  Table  3.4b. 

Table  3.4a 

Values  of  Coefficients  in 


Policy  p/h 

a 

2, 

0 /p 

b 

= 3 

c 

d 

u 

a 

2, 

0 /p 

b 

« 9 

c 

d 

Optimal  4 

.7164 

.5123 

. 3958 

.1614 

3.029 

.4356 

.4267 

-.0369 

9 

1.747 

.4236 

.3115 

-.0175 

3.550 

.3113 

.3072 

.0345 

24 

1.455 

.3525 

.3001 

.0563 

4.870 

.2512 

.2355 

.0067 

99 

2.031 

.2851 

.2660 

5.643 

. 1975 

.1986 

.0035 

Power  4 

.8282 

.3927 

.2867 

.2349 

.9887 

.4456 

.3352 

.2714 

Approx.  9 

.4450 

.4280 

.4255 

.3302 

1.722 

.4233 

.3892 

.1332 

24 

1.171 

. 3775 

. 3938 

.0932 

3.115 

.2864 

.3250 

.0814 

99 

1.015 

.2790 

4.237 

.0288 

-.0093 

.2013 
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Table  B.Ab 


Frequencies  of  A(Bj")  in  a 384-ltem  System 


Range 

for  A(B^') 

Number  of  Items 

Cumulative  Percentage 

I terns 
^wi  th 
o^/u=3 

Items 

-With 

o^/li=9 

A1 1 
Items 

1 1 ems 
-With 

0 /u=i 

of  Items 

Items 

2with 

0 /p-9 

All 

Items 

[0%.2%) 

34 

74 

108 

18% 

39% 

28% 

[2%. 4%) 

42 

55 

97 

40% 

67% 

53% 

[4%. 6%) 

38 

34 

72 

59% 

85% 

72% 

[6%. 8%) 

25 

14 

39 

72% 

92% 

82% 

[8%, 10%) 

29 

7 

36 

88% 

96% 

92% 

[l0%,12%) 

10 

6 

16 

93% 

99% 

96% 

[12%. 14%) 

5 

2 

7 

95% 

100% 

98% 

[14%, 16%) 

5 

5 

98% 

99% 

[l6%,18%) 

3 

3 

99% 

100% 

[l8%,20%) 

1 

1 

100% 

100% 

4 3h 

( 
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We  examined  other  ways  of  approximating  backlog  cost  with  multi- 
plicative expressions.  The  quantity  Il-B/(pp)]  was  fit  to  our  data 
without  success.  Even  with  separate  fits  for  each  value  of  p/h  our 


coefficients  of  determination  were  only  of  the  order  0.93.  We  also 
examined  the  quantity 

“f  "t  • 

where  , and  are  our  fixed-policy  multiplicative  expressions 

for  total  cost,  replenishment  cost,  and  holding  cost,  respectively.  This 
approach  also  was  unsuccessful;  deviates  from  B by  an  average  of  17%. 

3.2.3  Backlog  Protection 

To  obtain  a multiplicative  expression  for  backlog  protection,  we 
used  the  quantity  (l+p/h)P  as  the  dependent  variable  in  a regression 
model  of  the  form  (29).  (30).  (31).  For  optimal  policies  we  obtain  a 
coefficient  of  determination  of  0.99995.  and,  dividing  by  (1+p/h)  , 
we  have 

(3A)  - A9.A7  (p/h)“^’^^^^^Pexp(.010A5p/h)/(l+p/h)  . 


For  Power  Approximation  policies  we  have  a coefficient  of  determination 
of  0.99986  with  the  expression 


(35)  Pp  - A9.7A  (p/h)“^’^^^^'^Pexp(.010Alp/h)/(l+p/h)  . 

We  define  &(Pj)  as  the  percentage  deviation  of  (3A ) or  (35)  from 
the  actual  value  of  P for  an  item.  Table  3.5  is  a summary  of  the 
distribution  of  A(Pj.)  in  the  576-ltem  system.  The  average  value  of 

A(Pj)  in  the  system  is  0.7%  with  fewer  than  1%  of  the  items  deviating 


AA 


Table  3.5 


► 


;4»i 


Frequencies  of  A(P^)  in  a 576-item  System 


Range  for 

A(P^) 

Number 
of  Items 

Cumulative  Percentage 
of  Items 

[0%,2%) 

525 

91% 

(25!, 4%) 

43 

99% 

[4%, 6%) 

7 

ioo"% 

[6%, 85!) 

1 

100% 

[ 


I 


Vi 

V 

t 


by  more  than  4%.  The  pattern  of  errors  is  quite  similar  to  that  noted 

for  the  analytic  expressions  P in  Section  2.4.  The  largest  deviations 

a 

occur  for  items  with  small  p/h  . Items  with  p/h  equal  to  4 account  for 
all  errors  greater  than  4%,  and  44  of  51  items  with  deviations  greater 
than  2%. 

The  pattern  of  errors  for  is  nearly  identical  to  the  pattern 

noted  for  P (23),  (24)  in  Section  2.4  Since  P is  much  easier  to 
a a 

compute,  we  suggest  that  it  be  used  rather  than  . 

3.2.4  Average  Replenishment  Cost  per  Period 

The  following  approximation  for  optimal  replenishment  cost  yields  a 
coefficient  of  determination  of  0.996  and  an  average  deviation  of  2.25! 
from  actual  values  of  R 


R 

o 


.6050h  exp(-.08549a2/u)p‘ 


5126+.004333o2/p 


^L+i)--03847-.006147o2/p 

^j^^j^j.5107+.01492o2/p 


, ,,,,-.09095h/p+.02130p/o2 
(p/h) 


Table  3.6  is  a summary  of  the  distribution  of  percentage  errors  for  R^ 
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Table  3.6 


Frequencies  of  A(R  ) In  a 288-ltera  System 

Q 


Range  for 

A(R  ) 
u 

Number 
of  Items 



Cumulative  Percentage 
of  Items 

[0%,2%) 

163 

57% 

1 2%, 4%) 

88 

87% 

[4%, 6%) 

21 

94% 

[6%, 8%) 

11 

98% 

[8%, 10%) 

3 

99% 

110%, 12%) 

0 

99% 

[12%, 14%) 

1 

ioo‘% 

[14%, 16%) 

- 

1 

100% 

In  the  288-lteiti  system.  The  two  items  with  values  of  exceeding 

lOX  have  a variance  to-mean  ratio  9,  ieadtime  0,  mean  2,  and  unit  penalty 
cost  ratio  and  unit  setup  cost  ratio  of  4,64  and  99,32,  respectively. 

We  obtain  an  even  better  fit  for  Power  Approximation  replenishment 
cost.  The  follov^rg  approximation  has  a coefficient  of  determination  of 
0.999  and  an  average  percentage  deviation  of  1.4X  from  actual  values  of  R 

n ccou/  2/  V--2774  .4959+.007906a2/y-.003856p/o2 

R . 5763h(a'‘/p)  p 

P 

(L+l)--°57^^K/h)-^0^^-‘0^563p/a2  _ 

There  are  no  significant  outliers  in  the  distribution  of 

find  74%  of  the  288  Items  with  A(R  ) in  [0%,2X)  and  all  the  remaining 

P 

items  have  values  of  6(R  ) in  [2%, 4%)  . 

P 
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The  reason  for  R having  smaller  errors  than  R Is  tliat  the 
P o 

Power  Approximation  has  a simple  multiplicative  formula  for  setting 
(S-s).  The  multiplicative  nature  of  the  formula  and  its  Independence 
of  the  parameter  p/h  make  the  associated  cost  a natural  candidate  for 
a regression  fit  of  this  type. 

J . 2 . 5 Average  Total  Cost  per  Period 

We  would  expect  to  obtain  good  multiplicative  approximations  ♦'or 
total  cost  because  of  the  empirical  observations  described  in  Section  3.1 
We  also  would  expect  the  Power  Approximation  and  optimal  policy  fits  to  be 
similar  because  actual  total  cost  values  for  these  policies  are  nearly 
equal.  Both  of  these  conjectures  are  correct. 

The  following  approximation  for  optimal  total  cost  yields  a coeffi- 
cient of  determination  of  0.994  and  an  average  percentage  deviation  of  2.7% 

(36)  T^  = 2.663h  exp  ( . 097840^/^)  y ' 

(L^lj.l479+.01268a2/y-.06757y/o2 

^p^l^^-1.143h/p-.O2488y/02+.O13O2o2/y 

(j./h).4093-.02376o2/y 

The  following  approximation  for  Power  Approximation  total  cost  yields  a 
coefficient  of  determination  of  0.994  and  an  average  percentage  deviation 
of  2.8% 

(37)  9230h/p3-.0l508o*/w 

(L^lj.l498+.01231o2/y-.070j0y/a2 

. 3095-.01310o^/y+.  1073y/o^ 
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The  distributions  of  percentage  deviations  for  the  two  approximations 
are  nearly  identical.  We  summarize  them  as  om  system  in  Table  3.7  The 


Table  3.7 

Frequencies  of  A(T^)  in  a 576-item  System 


Range  for 
A(T^) 

Number 
of  Items 

Cumulative  Percentage 
of  Items 

[0%,2%) 

272 

47% 

[2%, 4%) 

176 

77% 

[4%, 6%) 

69 

90% 

(6%, 8%) 

31 

95% 

[SX.lOX) 

22 

99% 

[10%, 12%) 

1 

99% 

[12%, 14%) 

5 

100% 

largest  errors  are  consistently  associated  with  small  leadtime  and  large 
varlance-to-mean  ratio. 

We  attempted  to  Improve  tlie  approximation  by  fitting  separately  for 
each  value  of  o^/u  , but  the  approach  was  unsuccessful. 

The  approximation  does  Improve  when  we  fit  separately  for  each  value 
of  leadtime.  We  obtain  coefficients  of  variation  averaging  0.996  and  an 
average  percentage  deviation  of  2.1%  with  expressions  of  the  form 
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c +C  o^/u 

(38)  T^/h  = a expCbjO^/u+b^u/o^)^ 

d h/p+d  o^/u+cl,u/o^  e +e,o^/u 
(p/h)  ^ ^ ■’  (K/h)  ^ ^ 

Values  of  Che  coefficients  in  (38)  are  given  in  Table  3.7  for  each  policy 
and  each  value  of  L . We  summarize  the  distribution  of  errors  for  (38) 


Table  3.7 

Coefficient  Values  in  Expression  (61)  for  T^ 


Policy 

i 

IB 

*’2 

B 

"2 

B 

^2 

"l 

®2 

Optimal 

2.084 

.0798 

-.0153 

.4929 

-.0033 

-.8398 

.0158 

-.0198 

.4601 

-.0218 

I 

3.212 

.1196 

-.0584 

.4955 

-.0054 

-2.226 

.0124 

-.0282 

.4024 

-.0249 

1 

4.216 

.1289 

-.1042 

.4958 

-.0057 

-1.377 

.0108 

-.0280 

.3652 

-.0246 

2.095 

.0791 

-.0198 

.4951 

-.0041 

-.8052 

.0160 

-.0159 

.4552 

-.0212 

■ 

B 

3.216 

.1198 

-.0606 

.4976 

-.0058 

-1.211 

.0125 

-.0266 

.4003 

-.0248 

■ 

B 

4.216 

.1296 

-.1098 

.4968 

-.0058 

-1.358 

.0109 

-.0258 

.3644 

--  ■ 

-.0248 

in  Table  3.8.  The  number  of  items  having  deviations  within  2%  has  Increased 
by  more  than  a third,  and  the  number  of  items  with  deviations  exceeding  6% 
has  been  nearly  halved.  The  four  items  with  deviations  exceeding  lOX  are 
listed  In  Table  3.9. 

3. 3 Approximations  for  Statistical  Policies 

We  consider  four  72-ltem  systems  (see  Table  3.1)  controlled  with 
the  Statistical  Power  Approximation  (see  Section  1).  One  system  has  a t 
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Table  3.8 


Frequencies  of  A(Tj)  in  a 576-item  System 


Range  f or 

&(Tp 

Number 
of  Items 

Cumulative  Percentage 
of  Items 

10X.2X) 

364 

632 

122. A?) 

145 

882 

IA2.6X) 

36 

952 

(62.81) 

13 

972 

(H2.102) 

14 

992 

[102.1221 

4 

1002 

Table  3,9 

Items  with  A(Tp  Exceeding  10% 


Policy 

o^/u 

D 

D 

p/h 

K/h 

Error 

Optimal 

9 

0 

2 

99 

32 

-11.82 

Power  Approximation 

9 

0 

2 

99 

32 

-11.52 

Po\#er  Approximation 

9 

0 

2 

4 

32 

10.92 

Optimal 

9 

0 

2 

4 

32 

10.92 
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var lance-to-mean  ratio  of  3,  and  revision  interval  length  equal  to  26 
periods.  The  other  three  systems  all  have  a varlance-to-mean  ratio  of  9, 
with  revision  Interval  lengths  of  13,  26,  and  52  periods. 

We  develop  separate  operating  characteristic  approximations  for 
each  of  the  systems  using  regression  models  of  the  form  (29),  (30),  (32). 
We  label  a quantity  with  subscript  (x,y)  when  it  pertains  to  a system 
with  a^/p  equal  x , and  a y-perlod  revision  Interval.  For  example 
H _ is  the  holding-cost  approximation  for  the  system  with  a^/p  equal 

i f Zb 

3,  and  revision  Interval  equal  26  periods. 

We  also  develop  approximations  that  explicitly  include  revision 
interval  length  T as  a variable.  For  these  approximations  we  use 
regression  models  of  the  form  (29).  (30).  (33)  on  the  systems  with 
varlance-to-mean  equal  9.  We  label  these  approximations  with  the  sub- 
script (9,T)  . 

When  we  summarize  error  distributions  we  combine  the  systems  into 
one  288-ltem  system.  The  subscript  "s"  is  used  to  label  the  system. 

Recall  that  our  statistical  policy  data  contain  random  simulation 
errors.  Hence,  we  can  expect  that  the  quality  of  the  regression  fits 
will  degrade  as  compared  with  those  for  fixed  policies. 

3.3.1  Average  Holding  Cost  per  Period 

We  obtain  tiie  following  approximations  for  statistical  policy  holding 

cost 
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= 10.27hp-'^^^®(L+l)-^^^^(p/h)"^-^^^^''P(K/h)-°^^^^  . 

The  approximations  have  coefficients  of  determination  of  0.989,  0.990, 
0.992,  and  0.994,  respectively.  Average  percentage  deviations  from 
actual  viilues  of  H are  5.0%,  5.4%,  4.7%,  and  4.3%,  respectively,  as 
compared  with  3.9%  for  Power  Approximation  fixed  policies. 

The  percentage  deviations  are  distributed  similarly  for  the  four 
systems.  We  summarize  the  deviations  in  Table  3.10.  Unlike  the  fixed- 
policy  fits,  we  find  no  general  pattern  in  the  parameter  settings  of 

Table  3.10 

Frequencies  of  if'  ^ 288-ltem  System 


Range  for 

A(H  ) 
s 

Number 
of  Items 

Cumulative  Percentage 
of  Items 

[0%,2%) 

63 

22% 

(2%, 4%) 

76 

50% 

(4%, 6%) 

59 

66% 

[6%, 8%) 

43 

81% 

[8%, 10%) 

23 

89% 

[10%, 12%) 

11 

97% 

[12%, 14%) 

12 

ioo"% 

[14%, 16%) 

1 

100% 
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outliers.  Nearly  half  of  the  items  with  deviations  larger  tlian  10%  have 
p/h  equal  to  99. 


We  also  performed  separate  fits  for  each  value  of  p/h  , and  explicitly 
included  revision  interval  length  T as  a variable  in  the  o^/u  equal  9 
system.  The  average  percent  deviation  decreases  to  3.8%  for  the 
approximations 


H’ 


(3,261 


1.365hp' 

•207^K/h) 

.2319 

, p/h  = 4 

2.156bp 

•2^^\K/h) 

.2082 

, p/h  = 9 

A.792hp 

•3^««(K/h) 

. 1400 

, p/h  = 99 

Hi 


(9.T) 


2 . 7 /Ohp ' 

.4655+.00033T^j^^^^ 

.2902+1, 

.1155 

, p/h  = 4 

4 . OOOhp • 

,4587+.00042T^^^^^ , 

.2994+1. 

.1280 

, p/h  = 9 

9.698hp 

.4654+.00047T^^^^^ 

. 3935+1. 

.0566 

. p/h  = 99 

Although  the  stepwise  regression  program  permitted  the  variables  T and  1/T 
to  enter  every  exponent,  notice  tliat  T and  1/T  each  enter  only  once 
in  H’  . 

We  summarize  the  distribution  of  errors  for  the  approximations  in 
Table  3.11.  The  distribution  is  closer  to  0%  than  in  Table  3.10,  with 
27%  more  items  with  deviations  less  than  A%. 

3 . 3.2  Average  Backlog  Cost  per  Period 

We  encountered  serious  difficulties  In  fitting  multiplicative  expres- 
sions to  statistical  policy  backlog  costs.  When  we  use  the  model  (29), 
(30),  (32)  on  the  four  statistical  systems  we  obtain  coefficients  of 
determination  averaging  0.80. 
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Table  3.11 


Frequencies  of  A(ir)  in  a 288-item  System 


Range  for 

A(H’) 

s 

Number 
of  Items 

Cumulative  Percentage 
of  Items 

[0%,2%1 

96 

33% 

[2%, 4%) 

80 

61% 

[4%, 6%) 

48 

78% 

[6%, 8%) 

36 

90% 

[8%, 10%) 

L8 

97% 

[10%. 12%) 

5 

98% 

[12%, 14%) 

4 

ioo"% 

(14%,lf>%) 

1 

100% 
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We  obtain  improved  results  by  fitting  each  system  separately  for 
each  value  of  p/h  , yielding  expressions  of  the  form 

b -b  /p  c +c  (L+1)  dh/K 

= a h exp(bj^/M)p  (L+1)  (K/h) 

In  Table  3.12  we  list  values  of  the  constants  a,...,d  for  each 
statistical  system  and  each  value  of  p/h  . The  quality  of  the  fit  is 


Table  3.12 

Values  of  Constants  in  B 
s 


System  p/h 

a 

^2 

‘’3 

"i 

d 

(3,26)  4 

6,355 

-.6755 

.2258 

-1.067 

.3691 

.0087 

-3.989 

9 

2.599 

-.3709 

.3581 

.0952 

.4952 

.0066 

-2.697 

99 

.4842 

1.568 

,6172 

2.856 

.3656 

-.0132 

-1.380 

(9,13)  4 

1.983 

-.1610 

.5620 

1.012 

.6015 

-.0171 

-1.609 

9 

3,319 

-.3796 

.4052 

.4183 

.5499 

.0070 

-.3021 

99 

.1781 

2.833 

1.294 

6.461 

.3405 

-.0008 

-.6803 

(9,26)  4 

7.427 

,8536 

.2346 

-.7821 

.5433 

-.0105 

-2.791 

9 

1.380 

.4098 

.6079 

1.464 

.5070 

.0066 

-.2855 

99 

4.251 

.8241 

.2755 

2.513 

.1585 

.0207 

-1.120 

(9,52)  4 

4.428 

-.4380 

.3681 

-.1355 

.5166 

-.0126 

-2.763 

9 

4.376 

-.2043 

.3079 

-.0763 

.4799 

-.0031 

-1.333 

99 

7.242 

.5117 

. 06  30 

J 

1.690 

i 

.0555 

.0082 

-1.030 

y 


rather  good  for  p/h  equal  4 and  9 , with  an  average  deviation  of  4X 
from  actual  values  of  B . The  approximation  for  p/h  equal  99  , 
however,  deviates  by  an  average  of  10%.  We  summarize  the  distribution  of 
deviations  in  Table  3.13. 


Table  3.13 

Frequencies  of  A(B  ) In  a 288-item  System 


Range  for 
A(B^) 

Number  of  Items 

Cumulative 

Percentage  of 

1 terns 

Items 

with 

p/h<99 

Items 

with 

p/h=99 

All 

I terns 

I terns 

with 

p/h<99 

Items 

with 

p/h«99 

All 

1 1 ems 

10%.2%) 

72 

10 

82 

38% 

10% 

28% 

[2%, 4%) 

48 

11 

59 

62% 

22% 

49% 

14%, 6%) 

37 

8 

45 

82% 

30% 

65% 

[6%, 8%) 

16 

9 

25 

90% 

40% 

7 3% 

18%, 10%) 

6 

12 

18 

93% 

52% 

80% 

(10%, 20%) 

13 

36 

49 

100% 

90% 

97% 

[20%, 30%) 

9 

9 

99% 

ioo"% 

[30%, 40%) 

1 

1 

100% 

100% 

We  caution  that  the  approximations  may  not  be  robust.  One  reason 
Is  that  we  used  only  24  items  to  fit  the  seven  constants  in  each  regression. 
Another  reason  Is  that  most  of  the  coefficients  in  Table  3.12  do  not  vary 
monotonlcal ly  with  the  three  values  for  p/h  ; interpolating  coefficient 
values  for  other  settings  of  p/h  would  not  be  warranted. 
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3.3. 3 Backlog  Pr(  tectlon 


Recall  our  observation  In  Section  3.2.3  that  for  fixed-policy  backlog 
protection,  the  easy-to-calculate  expressions  (23)  and  (24)  are  as  accurate 
as  the  complicated  multiplicative  expressions  (34)  and  (35 j.  Therefore  we 
fit  functions  of  the  form  (23)  and  (24)  to  the  statistical-policy  data. 

We  have  average  percent  deviations  of  only  1.3%  for  the  following 
backlog  protection  approximations 

= (-.0925+.9995p/h)/(l+p/h) 

13, /o; 

1>(9  = (-.0057+.9853p/h)/(l+p/h) 

l’(9  263  = (.0l67+.9904p/h)/(l+p/h) 

H(9  32)  = (•0511+.9943p/h)/(l+p/h)  . 

We  summarize  the  distribution  of  errors  in  Table  3.14.  We  see  the  same 
pattern  of  errors  as  with  the  approximations  for  backlog  protection  in 
Sections  2.4  and  3.2.3.  The  largest  errors  are  due  to  Items  with  low 
values  of  p/h  , for  example,  fifteen  of  the  twenty  items  with  deviations 
greater  than  4%  have  p/h  equal  four  or  nine. 


Table  3.14 

Frequencies  of  A(P  ) in  a 288-ltem  System 


Range  for 

A(P  ) 
s 

Number 
of  Items 

Cumulative  Percentage 
of  Items 

10%,2%) 

223 

77% 

[2%, 4%) 

45 

93% 

14%, 6%) 

18 

99% 

[6%, 8%) 

2 

100% 
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The  approximation  Is  Improved  by  fitting  to  minimize  relative  squared 

error  as  in  Section  2.4.  Revision  interval  length  T is  also  Introduced 

2 

as  a variable  iti  the  systems  having  o /p  = 9.  The  result  is  the  following 
approximation  which  has  an  average  percentage  error  of  1.3Z 

P/o  = (-.272  + 5.66  h/K  + 1.003  p/h)/(l  + p/h) 

( d, 26) 

^(9  T)  ° + .994  p/h)/(l  + p/h)  . 

We  summarize  the  distribution  of  errors  for  this  approximation  in  T.ihle 
3.14a.  Although  1*'  has  the  same  average  percentage  error  in  the  28H-item 


Tab le  3.14a 

Frequencies  o£  ^(Pg^  ^ 288-ltem  System 


Range  for 

A(P') 

s 

Number 
of  Items 

Cumulative  Percentage 
of  Items 

[0%,2%) 

220 

76% 

[2%, 4%  ) 

53 

95% 

[4%,6X) 

15 

100% 

system,  there  are  fewer  items  with  large  percentage  errors  using 
approximation 
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3.3.4  Average  Replenishment  Cost  per  Period 

We  obtain  the  following  approximations  for  statistical  policy  replen- 
ishment cost 

R,  = .4313hp-'''\Lfl)--°''«^K/h)-^«18 

J , zo 

Rg  = .2676hp'^^®®(L+l)'''"’^®^(K/h)'^^^° 

R9  25  = .2949hy‘^^^^(L+l)"'°^^^^(K/h)'^^^^ 

R9  52  “ .3012hp  (L+1)  (K/h) 

The  approximations  have  coefficients  of  determination  of  0.999,  0.989, 

0.998,  and  0.999,  respectively.  Average  percentage  deviations  from 
actual  values  of  R ere  1.1%,  4.0%,  1.6%,  and  1.2%,  respectively,  as 
compared  with  2.0%  for  Power  Approxlmat ion  fixed  policies. 

We  summarize  the  distribution  of  errors  in  Table  3.15.  Of  the  34 
items  with  deviations  larger  than  4%,  28  are  from  the  (9,13)  system.  All 
of  the  items  with  deviations  larger  than  10%  are  from  this  system,  and  also 
have  L equal  4 , and  p/h  equal  99  . 

We  also  performed  separate  fits  for  each  value  of  p/h  , and  explicitly 
Included  revision  interval  length  T and  its  reciprocal  1/T  as  variables 
in  the  a^/y  equal  9 system.  The  average  percent  deviation  decreases 
from  2.0%  to  1.6%  for  the  approximations 


.4427hy'^^^^(L+l)" 

.0586 

(K/h) . 

p/h  « 4 

S 17  ^ 

.'.3641,, I fl.t  1; 

.051  . 

(l./h)-^^'^'^  . 

p/h  - 6 

. 1 '»4  , 

.0'i5  1 

. ..5603 

. 4 1 'j  ''yll.l  ( 1.  t 1 ) 

' 'll)  , 

P b ■ 

I 
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Table  3.15 


Frequencies  of  A(R  ) in  a 288-itera  System 


Range  for 

A(R  ) 
s 

Number 
of  Items 

Cumulative.  Percentage 
of  Items 

[03:,2%) 

202 

70% 

[2%. 4%) 

52 

88% 

[4%,  67.) 

20 

95% 

[6%, 8%) 

6 

97% 

[8X,102:) 

4 

99% 

[10%, 12%) 

1 

99% 

[12%, 14%) 

0 

99% 

[14%, 16%) 

2 

ioo"% 

[16%, 18%) 

0 

ioo"% 

[18%, 20%) 

0 

ioo"% 

[20%, 22%) 

0 

ioo“% 

[22%, 24%) 

0 

ioo"% 

[24%, 26%) 

1 

100% 

.2872hy'^^'*'°(L+l)”’°^^®~'^^^°^^(K/h)*^^^^  . p/h  - 9 

.2852hu*^^^^(L+l)''°^^^“^'^^^'^^(K/h) , p/h  - 99  . 

Observe  that  in  the  stepwise  least-squares  fitting  procedure,  the  revision 
Interval  entered  only  through  a single  term  In  1/T  , and  was  not  selected 
for  any  of  the  other  coefficients. 
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We  sununariEe  the  distribution  of  errors  for  the  approximations  in 
Table  J.16.  The  distribution  is  better  than  in  Table  d.l5,  with  outliers 
still  coming  from  the  (9,13)  system.  All  the  items  with  deviations  greater 
than  6%  come  from  this  system.  They  also  tend  to  have  L equal  4 , p/h 
equal  99  , and  small  u . 


r- 

1 


Table  3.16 


Frequencies  of  ^ 288-item  System 


Range  for 

a(r;) 

Number 
of  Items 

Cumulative  Percentage 
of  Items 

[0%,2%) 

208 

72% 

[2%, 4%) 

63 

94% 

[4%, 6%) 

11 

98% 

[b%,m 

2 

99% 

[8%, 10%) 

2 

99% 

[10%, 12%) 

1 

100~% 

[12%, 14%) 

0 

100~% 

[14%, 16%) 

1 

100% 

3.3.5  Average  Total 

Cost  per  Period 

We  obtain  the  following  approximations  for  statistical  policy  total  cost 


T 

T 

T 

T 


3.26 
9,13 

9.26 
9,52 


2.370hp 

■'''«(p/h) 

•'5^«(K/h) 

3.776hp' 

■““'(Wl)' 

•^^^Np/h) 

•''''(K/h) 

3.798hw 

■'■“’(LH) 

•'°2^p/h) 

•“'°(K/h) 

3.758hw 

•2^^^p/h) 

•^^^^(K/h) 

.3022 

.1843 

.1917 

.1964 
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The  approximations  have  coefficients  of  determination  of  0.989,  0.987, 


0.983,  and  0.985,  respectively.  Average  percentage  deviations  from 
actual  values  of  T are  3.8%,  5.1%,  5.4%,  and  4,9%,  respectively,  as 
compared  to  2.8%  for  the  Power  Approximation  fixed  policies. 

We  summarize  the  distribution  of  percentage  deviations  in  Table  3.17. 
The  patterns  of  percentage  errors  are  quite  similar  in  all  four  systems, 
with  large  errors  tending  to  occur  for  items  with  large  p/h  and  small  u • 

We  also  performed  separate  fits  for  each  value  of  p/h  , and 
explicitly  Included  revision  Interval  length  T as  a variable  in  the 
o^/u  equal  9 system.  The  average  percent  deviation  drops  from 


Table  3.17 

Frequencies  of  in  a 288-ltem  System 


Range  tor 

A(T  ) 
s 

Number 
of  items 

Cumulative  Percentagi 
of  Items 

10%,2%) 

83 

29% 

[2%, 4%) 

69 

1 

53% 

[4%, 6%) 

1 

39 

66% 

16%, 8%; 

43 

81% 

18%, 10%) 

19 

88% 

[10%, 12%) 

16 

93% 

[12%, 14%) 

11 

97% 

[14%. 16%) 

7 

ioo"% 

(16%, 18%) 

0 

ioo'% 

[18%, 20%) 

1 

100% 

4.8%  to  2.6%  for  the  following  approximations 


J 
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2.299h„-'’’^^(I.+l)-^^'‘''(K/h)-““  , p/h  - 4 

^b.26)  ■ , p/h  - 9 

6.607h„-''”^I,+l)-“''’(K/h)-^“^  , p/h  . 99 

3.790hp-‘’‘bL91)-““^-»‘”'bK/h)-“'"'““'°'f  . p/h  . 4 

t;,  ,,  - 6.002hp'‘'’“(L-hl)'^''”+-”“"(K/h)-'<’'"’^'‘“5''f  . p/h  - 9 

19.201..'^‘“a9l)-””«-ni/T„/„.0881*.4356/T  ^ „ 

We  summarize  the  distribution  of  errors  for  the  approximations  in 
Table  3.18.  The  distribution  is  markedly  better  than  that  in  Table  3.17. 

The  number  of  items  with  deviations  smaller  than  4%  has  Increased  by  mtsre 
than  half,  and  all  35  items  that  previously  had  deviations  greater  than  10% 
now  are  witliln  10%. 


Table  3.18 

Frequencies  of  ^ 288-item  System 


Range  for 

A(T’) 

s 

Number 
of  Items 

Cumulative  Percentage 
of  Items 

[0%,2%) 

130 

45% 

[2%. 4%) 

100 

80% 

(4%, 6%) 

49 

97% 

[6%, 8%) 

8 

ioo"% 

[8%, 10%) 

1 

100% 
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3^4  Measures  of  Robustness 

In  tills  section  we  analyze  the  robustness  of  our  multiplicative 
approximations  for  fixed-  and  statistical-policy  operating  characteristics. 
We  use  the  approach  of  Section  2.7  to  examine  the  aggregated  behavior  of 
the  approximations  and  to  investigate  their  accuracy  when  parameter  values 
arc  extrapolated. 

In  Tables  3.19  and  J.20  we  summarize  the  accuracy  of  our  fixed-policy 

approximations  11^  , , and  as  a function  of  the  parameter 

settings  in  the  576-item  system.  Each  entry  is  the  percentage  excess  of 

the  sum  of  approximate  values  over  the  sum  of  actual  values  for  all  items 

iiaving  the  given  parameter  setting.  As  we  noted  in  Section  2.7,  most  of 

the  large  errors  tend  to  cancel  when  the  items  are  aggregated  in  tiiis  way. 

Even  the  large  errors  we  found  for  backlog  cost  with  p/h  equal  99  do  not 

appear  in  Tables  3.19  and  3.20.  The  largest  systematic  errors  appear  for 

the  Power  Approximation  holding  cost  when  o^/y  equals  1 . This  pattern 

did  not  appear  in  Section  3.2  when  we  examined  single-item  behavior. 

Table  3.21a  Is  a multi-item  summary  of  errors  for  the  statistical- 

policy  approximations  ’ ^s  ’ ^s  ’ ^s  ’ ^s  ' summarize  the 

multi-item  behavior  of  H'  , R'  , and  T'  in  Table  3.21b.  These  tables 

s s s 

exhibit  the  same  property  noted  above  that  the  large  single-item  errors 

tend  to  cancel  when  the  items  are  aggregated.  Approximations  11^  , R^  , 

and  T'  appear  to  be  more  accurate  than  H , R , and  T , respectively, 
8 8 8 s 

in  all  systems  except  for  the  o^/p  equals  9 systems  where  R^  has 

■smaller  errors  than  R'  . The  largest  systematic  errors  appear  in  the 

s 

total  cost  approximations  T^  for  p/h  equal  to  99  and  p equal  to  2. 

We  test  the  robustness  of  our  fixed-policy  approximations  by  using 
them  in  the  64-ltem  system  of  Table  2.11.  The  parameter  settings  are 
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Percentage  Errors  of  Multiplicative  Approximations  for  a 288-ltem  System 
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Controlled  with  Statistical  Power  Approximation  Policies 
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Pfcrcentage  Errors  of  Multiplicative  Approximations  for  a 288-iteni  Systein 


7- 


Interpolated  and  extrapolated  in  such  a way  that  no  item  has  any  one 
parameter  with  the  same  value  as  In  the  576-ltem  system  used  to  derive 
the  approximations.  We  examine  errors  of  approximations  , and 

only.  We  already  have  a simple,  robust  approximation  for  P that  was 
analyzed  thoroughly  in  Section  2.7.  Our  approximation  for  B is  not 
tested  here  because  of  its  inaccuracies  and  lack  of  an  explicit  dependence 
on  the  parameter  p/h  . 

We  summarize  the  distributions  of  percent  deviations  of  , 

and  in  Table  3.22.  Optimal  and  Power  Approximation  policies  display 

similar  patterns  of  errors,  so  we  summarize  them  as  one  64-item  system. 

The  approximations  clearly  lack  the  robustness  that  we  noted  in  Section  2.7 
for  the  analytic  approximations.  Even  if  we  confine  our  attention  to  items 
with  fewer  than  two  extrapolated  parameters,  we  find  average  errors  of 
7.8%,  4.5%,  and  5.4%  for  holding  cost,  replenishment  cost,  and  total  cost. 

Table  3.22 

Percentage  Deviations  of  Multiplicative  Approximations 
in  a 64-item  System 


respectively.  The  worst  single-parameter  extrapolations  for  all  the 
approximations  are  0^/p  equal  to  15  and  u equal  to  0.5  . We  conclude 
that  the  multiplicative  approximations  should  be  used  only  with  great 
caution  when  parameters  are  extrapolated  beyond  the  ranges  of  the  576-item 
system  (Table  2.1). 

We  also  examine  the  behavior  of  the  statistical-policy  approximations 
Hs  • ^g  interpolated  and  extrapolated  parameter  values. 

We  are  unable  to  extrapolate  these  approximations  for  values  of  0^/p 
so  we  use  the  16-item  system  of  Table  2,11  with  o^/y  fixed  at  9 . The 
results  are  similar  to  the  fixed-policy  results.  In  Table  3.23  we  summarize 
the  distribution  of  errors  for  the  16-item  system  using  the  approximations 
for  a 26-pcriod  revision-interval  length.  Remember  that  the  deviations  are 
a composite  of  the  error  introduced  by  using  a multiplicative  form  as  well 

Table  3.23 

Percentage  Deviations  of  Statistical  Policy  Approximations 

in  a L6-itein  System 


as  a random  component  from  the  statistical  aspect  of  the  data.  Backlog 
protection  has  an  average  deviation  of  0.4Z.  If  we  consider  only  items 
with  fewer  than  two  extrapolated  parameters,  the  average  deviation  is 
only  0.22. 

The  holding  cost  approximation  has  an  average  deviation  of  13%, 

which  drops  to  4%  when  items  with  more  than  one  extrapolation  are  ignored. 
The  worst  single-parameter  extrapolation  is  for  p equal  0.5  . 

The  replenishment  cost  approximation  has  an  average  deviation  of 

18%.  If  we  consider  only  items  with  fewer  than  two  extrapolated  parameters, 
the  average  deviation  is  only  3%.  The  worst  deviations  occur  for  u 
equal  0.5  . 

I The  total  cost  approximation  T^  has  an  average  deviation  of  16%, 

\ which  drops  to  5%  when  only  single-parameter  extrapolations  are  included. 

The  worst  single-parameter  extrapolations  are  for  p equals  0.5  and  p/h 
equal  132. 

We  conclude  that  the  statistical-policy  approximations  should  not  be 
used  for  extremely  small  values  of  p , and  that  multiple-parameter 
extrapolations  should  be  considered  cautiously. 


f 
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4.  THE  SENSITIVITY  OF  OPERATING  CHARACTERISTICS 
TO  PARAMETER  VALUES 

In  SecLlons  2 and  3 we  developed  approximate  expressions  for  tlie 
operating  characteristics  of  fixed-policy  and  statistical-policy  in- 
ventory systems.  Now  that  we  have  the  approximations,  we  are  in  a 
position  to  perform  analyses  tliat  were  hitlierto  intractable.  In  tliis 
cha|)ter  we  will  make  simple  statements  about  the  following  questions: 

1.  Wliat  are  the  cost  and  service  implications  of  a change  in 
mean  demand? 

2.  What  are  the  effects  of  consolidating  demand  from  several 
distinct  warehouses  into  a single  central  warehouse? 

3.  By  liow  much  do  costs  rise  wlien  service  is  Increased? 

4.  How  much  Is  it  worth  to  obtain  quicker  delivery  of 
rt'pienishment  orders? 

5.  What  are  tiie  cost  implications  of  a change  In  unit  holding 
cost? 

6.  Wliat  are  the  effects  of  changing  the  review  period  length? 

For  each  question  we  present  calculations  based  on  a simplified 

model  to  Illustrate  the  uses  of  tiie  approximations  in  Section  3.  In 
many  cases,  our  assumptions  will  be  unrealistically  simple  for  the 
systems  planner  who  may  require  extensive  analyses  of  such  questions. 
Nevertheless,  our  results  indicate  general  trends  and  illustrate  the 
potential  uses  of  the  approximations. 

in  Sections  4.1  through  4.5  we  discuss  the  sensitivity  of  Power 
Approximation  and  Statistical  Power  Approximation  operating  character- 
istics to  mean  demand  p,  unit  penalty  cost  ratio  p/h,  leadtime  L,  setup 
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cost  ratio  K/h,  unit  holding  cost  h,  and  review  period  length.  In  each 
section  we  present  a simple  analysis  of  one  or  more  of  the  above  questions. 
We  also  include  graphs  displaying  the  operating  characteristics  of  Sec- 
tions 2and  3 as  functions  of  parameter  values.  Throughout  these  sections 
we  assume  that  o^/u  equals  9,  and  in  the  graphs  we  examine  single  para- 
meter variations  about  the  settings  m equals  8,  L equals  2,  p/h  equals 
9,  and  K/h  equals  A8. 

For  the  Power  Approximation  we  draw  graphs  using  the  analytic  ap- 
proximations (Section  2.3),  ^ (Section  2.4),  (Section  2.2), 

and  T (Section  2.5).  We  note  that  approximations  H , R , and  T are 

2 

discontinuous  functions  of  the  parameters  p,  o /m , L,  p,  K,  and  h.  This 

is  because  the  approximations  depend  on  the  policy  parameters  s and  U 

which  are  rounded  to  the  nearest  integer.  We  have  drawn  smooth  curves 

for  the  characteristics  by  omitting  the  round-off  procedure  for  s and  U 

in  the  computation  of  H , R , and  T . We  perform  sensitivity  analyses 

for  Power  Approximation  policies  using  approximations  (Section  3.2.1), 

P (Section  2.4),  R (Section  3.2.4),  and  T (Section  3.2.5).  With 
a,  p P P 

2 

o /m  equal  9 we  have 


(39) 


u 7 A7Q  K •'^629  ,,  ,s-3512  , -2.738  h/p  + . 0942 0691 

Hp  - 7.678  h M (L+  1)  (p/h)  (K/h) 

P - (.0695  + p/h)/(l  + p/h) 
a , p 

R - .3133  h (L  + l)"'°^^^(K/h)‘^^” 

P 

T 7 , -^^95  ,,  ^ ,,.2528  , ,.,-.9230  h/p  + . 1357 203 

Tp  - 5.663  h p (L  + 1)  (p/h)  (K/h) 


We  consider  the  statistical  policy  with  a 26-perlod  revision  interval 
2 

and  o /p  equal  9.  We  draw  graphs  using  two  sets  of  approximations.  When 

p/h  is  held  constant  at  9 we  use  H ' (Section  3.3.1),  R ' (Section  3.3.4), 

s s 
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(Section  3,3.5).  We  draw  the  remaining  graphs  and  perform 
sensitivity  analyses  using  the  following  approximations  from  Section  3.3 


A 


H 


(9,2b) 


10. AO  h u 


.4824 


(L  + 


(40) 


= ( -01^7  + .9904  p/h)  / (1  + p/h) 

(9, 2b) 

R = .2949  h (I.  + i)--0724  ^|^^,^^.62b9 

(9,2b) 


^ ,,  ...y  . .4309  ,,  ^ ,..  3024  , -,,.2550  ...,,x.l9l-  . 

T/o  = 3.798  h p (L  + 1)  (p/h)  (K/h) 

(9, 2bl 


We  do  not  use  Llie  approximations  for  backlog  cost  in  this  chapter 
because  of  tfieir  complexity  and  their  inaccuracy  for  certain  parameter 
settings. 


4 . 1 Sensitivity  to  Mean  Demand 

Consider  a situation  in  whicfi  mean  demand  cttanges  to  tfie  value 
(41  ) M2  = ^ 9i 

and  all  other  parameters  remain  unchanged.  From  (39)  and  (40)  we  see 
that  backlog  protection  will  be  approximately  unchanged  for  both  policies 
since  the  expressions  for  F are  independent  of  the  parameter  p.  All  other 
cost  expressions  are  of  the  form 

(42)  Cost  “ cp 

where  c and  u are  Independent  of  p.  I.et  Cost ^ and  Cost^  be  values  of  a 
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cost  component  when  p equals  pj  and  u 2 respectively.  Then  (41)  and  (42) 
yield 

Cost2/Costj  - a^, 

that  is,  the  ratio  of  approximate  costs  Is  a function  of  the  ratio  of  the 
two  mean  demands,  and  the  exponent  of  p in  the  cost  approximation.  In 
Table  4.1  we  display  approximate  percentage  changes  in  operating  characteris- 
tics for  several  values  of  the  parameter  a.  Similar  behavior  is  exhibited 
by  the  fixed  policy  and  the  statistical  policy.  All  three  cost  components 
are  roughly  proportional  to  the  square  root  of  p , with  replenishment  cost 

Table  4.1 

Approximate  Percentage  Changes  in 
Operating  Characteristics 


■ 

Policy 

Percentage 

Increase 

Holding 

Cost 

Backlog 

Protection 

Replenishment 

Cost 

Tota  1 

Cost 

in  M 

Power 

10% 

4.5% 

0% 

5.5% 

4.4% 

Approx. 

25% 

11% 

0% 

14% 

11% 

50% 

21% 

0% 

26% 

20% 

Statistical 

10% 

4.  7% 

0% 

5.4% 

4.1% 

Powe  r 

25% 

11% 

0% 

11% 

10% 

Approx. 

50% 

f 

22% 

0% 

25% 

1 

_ 

being  slightly  more  sensitive  than  the  other  costs  to  clianges  in  p . 

Now  we  consider  a situation  in  which  demand  from  n independent  and 
identical  warehouses  is  consolidated  into  a single  central  wareliousc. 
bet  Cost^  be  the  sum  of  values  for  a cost  i'ompt)nent  over  the  n wa r4'houst.'S 
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(before  consolidation),  and  Jet  Cost^  be  the  value  of  the  cost  component 
after  consolidation.  Using  (42)  we  have 

• u 

Cost  = n c p 

• (Jl 

Cost^  = c (n  p ) 

Combining  these  expressions,  we  have 

• f » “ L 

Cost2/CosLj^  = n 

which  is  independent  of  all  parameters  except  for  the  number  of  wan  - 
houses  before  consolidation,  and  the  exponent  of  p in  the  cost  approx ima L ion . 


Table  4.2 

Approximate  Percentage  Decreases  for 
Operating  Characteristics 


Policy 

.. 

■ 

Ho  Iding 

Backlog 

lep  lenishment 

Total 

Cost 

Protect  ion 

Cost 

Cost 

' 

i’owe  r 

2 

31% 

0% 

26% 

32% 

Approx. 

4 

52% 

0% 

45% 

53% 

8 

67% 

0% 

59% 

68% 

2 

30% 

0% 

27% 

33% 

Statist  lea  1 
Power 

4 

51% 

0% 

46% 

55% 

Approx. 

8 

J 

66% 

0% 

61% 

69% 
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Since  the  values  of  a in  (39)  and  (40)  are  all  less  than  1,  we  know  that 
Cost^  will  be  less  than  Cost^^  for  all  values  of  n greater  than  1.  In  n 
Table  4.2  we  list  approximate  percentage  decreases  in  costs  for  several 
values  of  n.  The  savings  are  considerable  for  both  fixed  and  statistical 
policies. 

Finally  we  consider  consolidating  demand  from  two  warehouses  that 
differ  in  one  respect.  Let  the  mean  demand  in  one  warehouse  be  mucli 
smaller  than  that  in  the  other  warehouse.  That  is,  the  item  is  a "slow 
mover"  in  one  location  and  not  in  the  otlier  location.  Let  the  ratio  of 
mean  demand  at  the  two  locations  be  a,  with  a less  than  one.  We  use 
(41)  to  calculate  the  total  cost  T^^  at  the  two  wareliouses  before 
consolidation,  yielding 

T^  = c[  + (a  M)“]  . 

After  consolidation  we  hove  total  cost  T^  given  by 

T2  = c [(1  + a)u]“  • 

After  combining  the  two  expressions  we  have 

T^/T^  = (1  + a) “/(I  + a“)  . 

We  also  consider  the  ratio  of  the  m.irginal  total  cost  of  the  sm.ill  ilenuind 
after  consolidation  to  its  total  cost  before  consolidation,  given 
approximately  by 

(T^-  CM'^*)/(a  m)““  [(!■  + a)'^-  ij/a** 

In  Table  4.3  we  list  approximate  values  for  total  cost  reductions  for 
several  values  of  a when  the  Power  Approximation  is  used.  The  reduction  in 
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Table  4.3 


Approxinvate  Percentage  Reductions  In 
Power  Approximation  ToLaJ  Cost 


Katio  of  Small  Demand 

Combined 

Marginal  Cost 

to  Large  Doraaiul 

System 

of  Slow  Movi'r 

0.1 

2 3% 

88% 

0,.3 

28% 

80% 

0.50 

31% 

73% 

total  costs  for  the  combined  system  is  quite  large  even  when  the  ratio  ol 
mean  denunds  is  only  0.1.  The  percentage  reduction  in  the  nvirginal  tost 
of  the  slow  iiwver  is  quite  draimitic,  especially  for  the  lowest  value  of 

I 

1 a.  The  results  for  statistical  policies  are  nearly  identical. 

in  Figures  4.1,  4.2,  and  4.3  we  present  graplis  of  the  approximations 
tor  holding  cost,  replenishment  cost,  and  total  cost,  respectively,  as 
functions  ot  mean  deiiund.  Tlie  fixed-  and  statistical-policy  a()proxi- 
mations  exhibit  similar  behavior.  The  statistical  policy  has  higher 
values  and  slightly  greater  sensitivity  to  p for  holding  cost  and  total 
cost.  Kep  I en  i shment  cost  values  are  nearly  Identical  for  both  pt)liiles. 
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roxlmaClons  versus  Mean  Demand 
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Replenishment  Cost  Approximations  versus  Mean  Demand 


MEAN  DEMAND 
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4.2  Sensitivity  to  Unit  Penalty  Cost 

Conslilvi  a situation  In  which  a systems  designer  wants  to  asscs.s  the 
cost  implications  of  Increasing  the  level  of  service.  The  designer  does 
this  by  assigning  a value  to  backlog  protection  P and  computing  the  corre- 
sponding value  of  the  unit  penalty  cost  ratio  p/h.  Approximate  values  lor 
the  operating  characteristics  are  then  calculated  using  this  value  of  p/h. 
Expressions  (39)  and  (40)  for  backlog  protection  are  of  tlu-  form 

P * (a  + b p/h)/(l  + p/h)  . 

We  solve  for  p/h,  yielding 

p/h  « (P  - a)/(b  - P)  . 

Suppose  backlog  protection  is  increased  from  P^  equals  0.90  to  P^  equals 
0.95.  For  the  I’ower  Approximation  We  have 

Pj^/h  =■  8.2 

p^/b  - 17.3, 

and  for  the  Statistical  Power  Approximation  we  have 

P^/h  i 9.8 
P^/h  - 23.1  . 

For  both  policies  the  unit  penalty  cost  ratio  increases  more  than  two-fold. 

We  list  tfie  approximate  percentage  Increases  in  the  operating 
characteristics  in  Table  4.4. 
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'I'jihli'  A.  A 


Approximate  Percentage  Increases  In  Operating  Characteristics 
When  P Increases  from  0.90  to  0.95 


Policy 

Holding 

Cost 

Backlog 

Protection 

leplenishment 

Cost 

Total 

Cost 

Power  Approximation 

38% 

6% 

0% 

20% 

i Statistical  Power 
, Approximation 

AA% 



6% 

0% 

2 5% 

Replenishment  costs  remain  approximately  unchanged  since  the  expressions 
for  this  characteristic  are  independent  of  p/h.  All  other  costs  rise  by 
20%  to  AA%,  with  the  statistical  policy  displaying  larger  increases. 

In  Figures  A.  A,  A.  5,  and  A.  6 we  present  graphs  of  the  a|)prox Inuit Ions  for 
lioldlng  cost,  backlog  protection,  and  total  cost,  respectively,  as 
functions  of  p/h.  All  the  curves  are  concave  over  the  range  of  values 
plotted.  The  holding  cost  curves  display  higher  values  for  the  statis- 
tical policy  than  for  the  fixed  policy.  The  curves  for  backlog  protection 
are  nearly  parallel,  with  larger  values  of  P for  the  fixed  policy.  For  both 
policies  we  see  that  backlog  protection  is  much  more  sensitive  to  p/li  for 
small  values  of  p/h.  The  total  cost  curves  show  that  the  statistical 
pollcv  has  higher  values  of  total  cost  and  greater  sensitivity  of  total 
cost  to  p/h  for  all  values  displayed  in  the  figure. 
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4 . 3 Sensitivity  to  Leadtime  and  Setup  Cost 


In  this  section  we  consider  the  fourth  question  raised  at  tlu-  begin- 
ning of  the  chapter.  Suppose  a quicker  means  of  replenishment  delivery  is 
being  considered.  What  is  it  worth?  We  assume  that  a quicker  delivery  is 
paid  for  by  means  of  a higher  setup  cost.  Therefore  we  are  confronted  with 
a trade-off  decision  where  leadtime  L is  decreased  and  setup  cost  K is 
increased . 

For  simplicity,  we  assume  that  the  time  between  ordering  and  deliv- 
ery (L^  + 1)  is  .in  even  number  of  periods  and  that  it  can  be  reduced  to 
(L^  + l)/2  periods.  The  setup  cost  is  simultaneously  Increased  from  Kj 
to  aKj^,  where  a is  greater  than  1.  We  assume  that  all  other  parameters 
remain  unchanged.  From  (39)  and  (40)  we  see  that  total  cost  per  period 
is  of  the  form 

T - c (1,  + 1)“  ( K/h)*^ 

where  c,  a , and  b are  independent  of  L and  K.  The  total  cost  Tj  lief  ore 
L and  K are  changed  is  approximated  by 

T^  - c (1.J  + O'*  (Kj/h)^^. 
and  after  L and  K are  changed  we  have 

T2  - c [(L^+  l)/2]“(aK^/h)^. 
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We  estimate  that  the  reduced  leadtime  will  realize  a savings  in  total 
cost  If 

T^/T^  = a^2"“  < 1 . 


In  other  words,  approximate  total  cost  is  reduced  if 


* 

a < a 


We  interpret  a as  an  approximate  break-even  point.  That  is,  if  the 

it 

setup  cost  corresponding  to  the  quicker  delivery  is  less  than  a K,  then 

the  quick  delivery  method  has  lower  total  cost.  Using  values  of  c»  and  li 

* 

from  (39)  and  (AO)  we  see  that  a equals  2. A for  the  Power  Approximation 
policy  and  3.0  for  the  statistical  policy.  Therefore,  it  is  worthwhile 
to  pay  more  t.han  double  the  setup  cost  for  a delivery  scheme  tliat  is 
twice  as  fast. 

We  calculate  numerical  results  for  the  case  of  a equals  2.  The 
cost  components  in  (39),  (AO)  are  all  of  the  form 

Cost  “ c (L  + 1 ) "^  (K/h)^  . 


Let  Costj^  be  the  value  of  a cost  component  before  changing  L and  K.,  and 
let  Cost^  i>e  the  value  after  the  change.  Then  we  have 

* B"*  Cl 

Cost2/Costj^  •=  2 


Specifically,  for  the  Power  Approximation  we  Iwve 


"2'"i 


.,-.282 


R2/RJ  « 2'^^^ 


• - 0A9 

T /T  - 2 ^ 

2'  1 
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and  for  the  statistical  policy  we  have 


H2/H^  » 2 


-.272 


R2/Rj  “ 2 


.699 


-.111 


» 2 


We  list  numerical  results  in  Table  4.5.  The  quicker  delivery  reduces 


Table  4.5 

Approximate  Percent  Changes  in  Operating  Characteristics 
L (L2  + 1)/(L^  + 1)  » 1/2;  K2/K^  - 2 


Policy 

Holding 

Cost 

Backlog 

Protection 

Replenishment 

Cost 

Total 

Cost 

Power  Approximation 

-18% 

0% 

57% 

-3% 

Statistical  Power  Approx. 

-17% 

0% 

62% 

-7% 

total  cost  by  3%  for  the  Power  Approximation  and  by  7%  for  the  statistical 
policy.  These  are  impressive  savings  since  we  assumed  that  the  quick 
delivery  was  twice  as  expensive  as  the  normal  delivery.  Backlog  protec- 
tion is  approximately  unchanged  since  the  expressions  for  P are  independent 
of  L and  K.  The  quick-delivery  system  has  lower  holding  costs  and  higher 
replenishment  costs. 

In  Figures  4.7,  4.8,  and  4.9  we  present  graphs  of  the  approximations 
for  holding  cost,  replenishment  cost,  and  total  cost,  respectively,  as 
functions  of  leadtime.  Although  the  curves  are  plotted  as  if  leadtlme 
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Power  Approximation 
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roxlmations  versus  Lead time 
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were  a continuous  variable,  we  note  that  leadtime  is  required  to  i)e 
integer-valued  in  our  model.  The  curves  for  holding  cost  and  replenish- 
ment cost  differ  only  slightly  for  fixed-  and  statistical-policy 
approximations.  Replenishment  cost  is  especially  insensitive  to  changes 
in  leadtime.  The  curves  for  total  cost  display  a greater  sensitivity 
to  leadtime  lor  the  statistical  policy  than  for  the  fixed  policy.  Tlie 
statistical  policy  curve  also  lias  larger  values  of  total  cost  for  all 
values  of  leadtime  except  for  a small  region  near  L equals  0.  In  tliis 
region  the  values  of  total  cost  may  be  considered  essentially  Lie.  ;;.nne 
since  the  difference  never  exceeds  2Z  which  is  less  than  the  average 
errors  of  approximation  for  total  cost. 

In  Figures  4.10,  4.11,  and  4.12  we  present  graphs  of  the  approxi- 
mations for  tiolding  cost,  replenishment  cost,  and  total  cost,  respectively, 
as  functions  of  K.  The  replenishment  cost  curves  are  nearly  idenli('al 
for  statistical-  and  fixed-policy  approximations.  For  liolding  cost,  Ltie 
statistical-policy  curve  has  higher  values  and  a steeper  slope  than  the 
fixed-policy  curve.  The  total  cost  curves  are  essentially  two  parallel 
lines,  with  the  statistical-policy  curve  about  four  units  above  the 
fixed-policy  curve. 
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Replenishment  Cost  Approximations  versus  Setup  Cost  Ratio 
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Total  Cost  Approximations  versus  Setup  Cost  Ratio 


r 3 


oo-gtT 
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4.4  Sensitivity  to  Unit  Holding  Cost 

Suppose  the  unit  holding  cost  changes  from  h^  to  h^.  What 
happens  to  the  operating  characteristics?  We  assume  that  the  de- 
signer chooses  to  maintain  the  same  Level  of  backlog  protection,  so 
p/ll  will  remain  unchanged.  From  (39)  and  (40)  we  see  that  all  tin.- 
other  characteristics  are  of  the  form 

Cost  = c h (K/h) ^ , 

where  c and  a are  independent  of  h.  Let  Cost^  and  Cost2  be  the  values 

of  a cost  component  when  h equals  h^  and  h2,  respectively.  We  then 

have  ^ ^ 

Cost2/Costj^  = (h2/hj^)  . 

Specifically,  for  the  Power  Approxinvit ion  we  have 

• .931 

H2/Hj^  = (h2/h^) 

R2/R^  = (h2/h^)‘^®^ 

T2/T^  » (h2/h^)-^^^, 

and  for  the  statistical  policies  we  have 
H2/H^  = (h2/h^)-^°° 

R2/R^  - (h2/h^)'^^^ 

T2/TJ  - (h2/h^)‘^‘^‘*. 

We  present  numerical  results  in  Table  4.6  for  the  case  of  h2/li| 
equal  1.1. 


96 


Table  4.6 


Approximate  Percentage  Increases  In  Costs 
(h^/h^  * 1.1;  p/h  held  constant) 


, - - 

1 

1 Policy 

i... 

Holding 

Replenishment  j 

Total 

Cost 

Cost  1 

Cost 

T 1 

Power 

9% 

4% 

8% 

Approximation 

, 

Statistical 

Power 

Approximation 

9% 

1 

4%  1 

i 

1 

8% 

In  Figures  4.13,  4.14,  and  4.15  we  present  graphs  of  the  approxl- 
nuilions  for  holding  cost,  rep lenisliment  cost,  and  total  cost,  respectively, 
as  functions  of  h.  As  in  the  discussion  above,  we  have  held  p/l>  constant 
at  9.  The  curves  for  holding  cost  and  replenishment  cost  are  nearly  the 
same  for  fixed  and  statistical  policies.  The  statistical-policy  curve 
for  total  cost  is  about  10%  higher  than  tlie  fixed-policy  curve.  Noli? 
that  the  curves  for  liolding  cost  and  total  cost  are  nearly  straight  lines. 
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4 . S Sens  it  to  Review  Period  Len^tli 


In  this  section  we  examine  the  efiects  of  changing  tiie  ri-vi<’w 
|>erio(J  ieiigth  liy  comparing  two  systems.  Let  tlie  review  period  in 
System  2 he  twice  as  long  as  tliat  in  System  1.  We  list  the  p.irameters 
ot  the  systi'iiis  in  'i'ahle  4,7. 


Tal)le  4.7 
Parameter  Values 


Sy  s ( iMii 

u 

2 , 

0 /|i 

p/ll 

(L  + 1) 

K 

h 

1 

9 

Pi/hf 

s 

s 

2 

9 

Pi /hi 

:Lj+1)/2 

*^1 

ah, 

The  denwiiui  per  period  in  System  2 is  a convolution  of  two  periods' 

2 

demand  in  System  1,  so  p is  twice  as  large  in  System  2 and  o /ii  Is  tlie 
same-  in  both  systems.  We  assume  that  the  systems  are  designi’d  for  the 
sam(>  level  of  backlog  protection  so  p/h  is  held  constant.  The  time 
between  placing  a replenishment  order  and  its  delivery  is  constant,  so 
we  set  it  at  D periods  in  System  1 and  (hj^+  •l)/2  periods  in 

System  2,  witli  1)  assumed  to  be  an  even  Integer.  The  setup  cost 

per  replenishment  order  K is  the  same  in  both  systems. 

We  assume  that  the  unit  holding  cost  hj  changes  by  a multiplicative 
constant  a when  the  review  period  lengtli  is  doubled.  We  expect  a to  take 
on  a value  between  1 and  2 in  most  systems.  It  would  equal  1 if  all  the 
holding  cost  were  attributable  to  review  functions  that  occur  once  each 


lOl 


period.  On  the  other  hand,  a would  equal  2 if  all  the  holding  cost  were 
proportional  to  time- in-storage. 


I 


We  now  examine  the  effect.s  of  review  period  length  on  the  cost 

2 

components  of  the  Power  Approximation.  Since  o /\i  and  p/h  are  the 
same  in  both  systems,  we  express  the  cost  approximations  (39),  (40) 
for  Systems  1 and  2 in  the  form 

Costj  - c h^M  " (L^+  1)*^  (K^/h^)* 

Cost^  » c(ah^)  (2m^)“[(L^+  D/d  ^ (K^/ahj)^ 

where  c, a , B , and  6 are  independent  of  review  period  length.  We 
take  the  ratio  of  costs  per  unit  time  by  dividing  the  costs  in  System  2 
by  twice  the  costs  in  System  1,  yielding 

Cost2/2  Cost^  = . 

Specifically,  we  have 

• 931. 

H2/2H^  = .540  a 

R2/2R^  = .771 

= .573  . 

In  Table  4.8  we  list  approximate  values  for  the  percentage  by  which  costs 
in  System  2 differ  from  those  in  System  l,for  several  values  of  a.  Sub- 
stantial cost  reductions  occur  for  the  longer  revision  Interval  when  a 
Is  close  to  1,  and  very  small  changes  occur  when  a equals  2. 

For  the  statistical  policies  we  assume  that  the  policy  ii'vlslt)n  In- 
terval is  the  same  length  of  time  in  both  systems.  That  is,  revision 
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Table  4.8 


Approximate  Percentage  Changes  In  Costs  per  lln i ii^ 
(Power  Approximation  Policies) 


.1 

Holding 

Cost 

Replenishment 

Cost 

Tota  1 

Cost 

1 

-46% 

-23% 

-43% 

1 . 5 

-21% 

- 9% 

-21% 

2.U 

3.0% 

1.9% 

5% 

occurs  every  T|  periods  in  System  1,  and  every  'i^/1  periods  iti  System  2, 
with  Tj  assumed  to  be  an  even  integer.  We  use  the  approximations  ob- 
tained in  Cliapter  9 tiiat  are  explicit  functions  of  revision  interval 
length.  We  si't  p/h  equal  to  9 and  use  the  approximations 


(41) 

'(9,T) 


4.000  h o-”"'  ♦ (I.  .4  I)-”’*  + ‘■““''‘(K/l,)-'-'™ 

.2872  1.  (I,  4 

6.002  I, <1.  4 ,)••*'’«  + “(K/l,)-™'’"  ■ 


We  express  these  cost  components  in  the  form 


Cost 


c 


li  P 


ori-rT 


(L  + 1) 


b+m/T 


(K/h) 


6+n/T 


where  c.oi,  b,  <5,r  , m,  and  n are  independent  of  revision  interval  length. 
As  before,  we  take  the  ratio  of  cost  in  System  2 to  twice  the  cost  in  Sys- 
tem 1,  yielding 


Cost^/^  Costj 


2 


o-b- 1 

3 


l-6j.  1^2 *■'^1''^  “ ^'"''’^1 

(Lj+  1) 


-2n/T, 
a I 

(Kj/hj) 


i 


Specifically,  we  have 


= 


R2/2K^ 


T2/2TJ  = 


.5584  a 


872 


7636  a 


.5752  a 


.367 


.00021T  -2.50/T  -.00021T,  ,,  ^ ,,1.250/T 

[2  1 1 1 (L^+  1) 


.799 


^-.223/T,  ^-..95/Tj  ^,.,76/r, 


We  evaluati-  a numerical  example  for  the  case  of  M ^ equals  4,  l.j  equals  3 
K^/h^  equals  64,  and  equals  26.  The  results  are  listed  In  Table  4.9 

for  three  values  of  a.  Each  entry  in  the  table  Is  the  percentage  by 


Table  4.9 

Approximate  Percentage  Changes  In  Cost  per  Unit  Time 
(Statistical  policy  with  M = 4,  » 3,  » 64,  Tj^  = 26) 


a 

Holding  1 
Cost 

Replenishment 

Cost 

Total 

Cost 

1 

-44% 

-24% 

-41% 

1.5 

-21% 

-11% 

-19% 

2.0 

1.8% 

-1.5% 

1.3% 

which  a cost  in  System  2 differs  from  a cost  in  System  1.  The  effeiits 
of  the  change  in  revision  interval  lengtli  are  nearly  the  .same  as  for 
the  Power  Approximation.  Substantial  cost  reductions  occur  lor  t he 
longer  revision  interval  when  a is  near  1,  and  very  small  changes  ociur 
when  a equals  2. 


I 12/T 
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4.6  Sensitivity  to  Variance-to-meaii  Ratio 

In  Figures  4.16,  4.17,  and  4.18  we  present  graphs  of  the  apf)roxl- 

iiiations  for  holding  cost,  replenishment  cost,  and  total  cost,  respectively, 

,is  functions  of  o /p  . We  have  curves  only  for  Power  Approximation  policies 

because  the  statistical-policy  approxiiuii t ions  do  not  have  an  ex|)lii  il  di'- 
2 

pendence  on  a /p . We  do  not  have  a curve  for  backlog  protection  since  t tu- 

approximation  for  P is  a function  of  only  p/h. 

The  curve  for  holding  cost  is  nearly  linear  with  an  increase  of  approx- 
2 

imately  70%  when  o /p  is  increased  from  one  to  nine.  The  repl  c-n  i shmen  t 

cost  curve  is  convex  in  shape,  decreasing  by  about  20%  when  o^/p  increases 

from  one  to  nine.  For  total  cost  we  see  a nearly  linear  relationship  and 

2 , 

an  increase  ol  about  50%  when  o /pis  increased  from  one  to  nine. 
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